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ABSTRACT 

We present the broad-band UV through mid-infrared spectral energy distributions (SEDs) of a sam¬ 
ple of 72 spectroscopically-confirmed star-forming galaxies at z = 2.30 ± 0.3. Located in a 72 arcmin 2 
field centered on the bright background QSO, HS1700+643, these galaxies were pre-selected to he at 
z ~ 2 based solely on their rest-frame UV colors and luminosities, and should be representative of 
UV-selected samples at high redshift. In addition to deep ground-based photometry spanning from 
0.35 — 2.15 fi m, we make use of Spitzer /IRAC data, which probes the rest-frame near-IR at z ~ 2. 
The range of stellar populations present in the sample is investigated with simple, single-component 
stellar population synthesis models. The inability to constrain the form of the star-formation history 
limits our ability to determine the parameters of extinction, age, and star-formation rate without 
using external multi-wavelength information. Emphasizing stellar mass estimates, which are much 
less affected by these uncertainties, we find (log = 10.32 ± 0.51 for the sample. The addition 
of Spitzer /IRAC data as a long-wavelength baseline reduces stellar mass uncertainties by a factor of 
1.5 — 2 relative to estimates based on optical - K s photometry alone. However, the total stellar mass 
estimated for the sample is remarkably insensitive to the inclusion of IRAC data. We find correlations 
between stellar mass and rest-frame R band (observed K s ) and rest-frame 1.4 /im (observed 4.5 fim) 
luminosities, though with significant scatter. Even at rest-frame 1.4/zm, the M/L ratio varies by a 
factor of 15 indicating that even the rest-frame near-IR, when taken alone, is a poor indicator of 
stellar mass in star-forming galaxies at z ~ 2. Allowing for the possibility of episodic star formation, 
we find that typical galaxies in our sample could contain up to three times more stellar mass in an old 
underlying burst than what was inferred from single-component modeling. In contrast, mass estimates 
for the most massive galaxies in the sample (M* > 1O 11 M 0 ) are fairly insensitive to the method used 
to model the stellar population. Galaxies in this massive tail, which are also the oldest objects in the 
sample, could plausibly evolve into the passive galaxies discovered at z ~ 1.7 with near-IR selection 
techniques. In the general framework of hierarchical galaxy formation and mergers, which implies 
episodic star-formation histories, galaxies at high redshift may pass in and out of UV-selected and 
near-IR color-selected samples as they evolve from phases of active star formation to quiescence and 
back again. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — galaxies: 


Within the last few years, new techniques and fa¬ 
cilities have resulted in a burst of activity focused on 
understanding the progress of galaxy formation dur¬ 
ing the era near z ~ 2. While this epoch has been 
kn own for som e time as the peak of the “quasar era” 
iSchmidt et all (1995)), until recently it has comprised 
part of the “redshift desert” about which little was known 
for normal galaxies. The recent exploration of galaxies 


starburst — stars: formation 
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in the z ~ 2 era has occurred on a number of fronts: 
new surveys have now successfully identified samples of 
2 ~ 2 galaxies based u pon t heir very bright far-IR and 
radio emission (e.g, iChaoma n et al .11 20041). their signa¬ 
ture rest-frame far -UV colors llAdelberger et alJl2Q04bt 
iSteidel e t al 120( 141) . their hig h luminosity in the observed 
K b and (|Daddi_el_a]j|2 004b) , red rest-frame op tical col¬ 
ors llvan Doldumretall2003t iFranxet aflEof) If), or some 
combination of the above (e.g.. lAbrahainet, alJ 12004 

iShaolev et al .1 l2004lf . The results of most of these sur¬ 
veys suggest that it was near 2 ~ 2 that galaxies as 
massive as any we see in the local universe began reach¬ 
ing maturity, and where a transition from rapidly star¬ 
forming to asymptotically quiescent may be occurring 
“before our eyes” for relatively massive galaxies (e.g., 
McCarthy et al.ll2004 IDaddi et al.H2004at IShaolev et al.l 

2004ft . Fortuitously, galaxies in the redshift range 
2 Sz z lij 2.5 are particularly well-suit ed to spectros copic 
studies both in the observed optical (ISteidel et al.ll2004t 
IDaddi et alJl2004bT) and especially in the observed near- 
IR, where nebular lines valuable for measurement of kine¬ 
matics, star-formation rates, and chemical abundances 
are observationally access i ble (e.g.j Erb et al.l2003ll2 004 
120051 IShaolev et alJl2004 Ivan Dokkum et al .112004) . 

One of the fundamental physical measurements that 
can be made for high redshift galaxies, and one of the 
more important for placing them into the theoretical con¬ 
text, is the stellar mass. A number of authors in the 
last several years have interpreted observed-frame UV 
to near-IR photometry with stellar population synthesis 
models in order to infer the stellar mass for samples of 
2 > 2 galaxies (e.g ..lPa povich et alJl2001t IShaolev et al.l 
120011: IDaddi et, all l2004bft . Others have used similar 
methods and the technique of photometric redshifts to 
obtain measu res of the total stellar mass as a function of 
redshift (e.g. iDickinson et al .1120014 iFontana et al .1120001: 

IClazebrook et al . iHf - Stellar mass measures become 

increasingly uncertain as one goes to higher redshift, 
where ground-based observations, limited to A 2.3 ^m, 
probe increasingly shorter rest-frame wavelengths. At 
2 > 2, the K s band corresponds to the rest-frame optical, 
where current or recent star formation may dominate the 
observed luminosity rather than the stellar mass that has 
accumulated over a galaxy’s lifetime. For this reason, the 
arriva l of data from the IRAC instrument llFazio et all 
l200l on board the Spitzer Space Telescope (hereafter 
Spitzer ) has been widely anticipated as a means to sig¬ 
nificantly improve the measurement of stellar mass and 
star formation history since it is capable of obtaining very 
sensitive observations out to 8/rm, allowing photometric 
access to the rest-frame 1 — 2pm light from galaxies to 
redshifts 2 ~ 4, and providing particularly sensitive ob¬ 
servations of this region for galaxies with more modest 
redshifts near 2 ~ 2. 

In this paper, we discuss the broad-band properties 
of a reasonably large sample of spectroscopically con¬ 
firmed, UV-selected galaxies at 2 ~ 2 in a field that was 
observed by the IRAC “In-Orbit Checkout” (IOC) pro¬ 
gram in the early days of Spitzer. The principal aim 
has been to explore the broad-band spectral energy dis¬ 
tributions (SEDs) for a robustly defined sample of high 
redshift galaxies that naturally spans a large range in 
stellar mass, star-formation rate, and star-formation his¬ 


tory. The main advantages of a UV-s elected sample at 
2 ~ 2, described in detail elsewhere iAdelberger et all 
12004b!: ISteidel et aiJl200l are that confirming optical 
spectra can be easily obtained, photometric selection is 
highly efficient and requires using only relatively “inex¬ 
pensive” optical ground-based photometry, and it is pos¬ 
sible to tune color criteria to yield a roughly volume- 
limited sample over the desired redshift range. In this 
particular case, an additional advantage was that the 
spectroscopic sample of galaxies already existed, as did 
deep photometry from observed 0.35 to 2.15 pm and tar¬ 
geted near-IR spectroscopic observations. Details of the 
ne ar-IR phot ometry and spectroscopy will be presented 
in i lErb et all (2005). The HS1700+6416 field is one of 
several we have used in a survey designed to combine ob¬ 
servations of 2 ~ 2 galaxies with sensitive observations 
of the intergalactic medium (IGM) in the same volume. 
In this respect, HS1700+6416 ( z em = 2.72) is one of the 
very best-studied high redshift QSOs in the entire sky, 
due to its unusual brightness (V ~ 16). It also happens 
to lie in the “continuous viewing zone” for both Spitzer 
and HST; this fact, plus the availability of the ancillary 
deep ground based data, is the main reason the field was 
chosen for the deep IOC observations by the IRAC team. 

In combination with the ground-based photometric 
and spectroscopic data, the IRAC observations of this 
field provide a means to establish the range of proper¬ 
ties represented among galaxies selected using UV colors; 
these may be useful for comparison to galaxy samples 
at similar redshifts but selected using other techniques. 
Of particular interest to us is a quantitative assessment 
of the “value added” by longer wavelength photometry 
for assessing the galaxy properties at high redshifts as 
compared to what one would infer from ground-based 
measurements alone (since realistically only a fraction of 
ground-based fields will be observed by Spitzer during its 
lifetime). As it turns out, the HS1700 field is also inter¬ 
esting because it allows for a preliminary comparison of 
the dependen ce of 2 ~ 2 galaxy p roperties on large-scale 
environment (ISteidel et a,1.11200511 . 

The paper is organized as follows: m presents the 
ground-based and Spitzer observations; describes the 
photometry and the selection of the sample which is dis¬ 
cussed in the rest of the paper; 0 describes the use 
of stellar population synthesis models to infer physical 
properties of the galaxies from their broad-band spectral 
energy distributions (SEDs); (O summarizes the results 
of the modeling focusing on the galaxy stellar masses 
and mass-to-light ratios (M/L); 0 discusses the results 
and their general implications; 32 summarizes the main 
results. 

We assume a cosmology with fl m = 0.3, Da = 0.7, and 
h = 0.7 throughout. 

2. OBSERVATIONS 
2.1. Optical Imaging 

Images of the HS1700+64 field were obtained in 2001 
May using the William Herschel 4.2m telescope on La 
Palma and the Prime Focus Imager. The WHT system 
uses a 2-chip mosaic of Marconi 2k x 4k CCDs with espe¬ 
cially good UV quantum efficiency, and provides a field 
of view of 16' by 16' sampled at 07236 per pixel. The 
conditions during the observing run were excellent, with 
sub-arcsecond seeing in all bands. The field was roughly 
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centered on the position of the QSO, and 15 arcsecond 
dithers were executed between subsequent exposures in 
order to allow for fringe removal in the 72 band and the 
creation of dark sky superflats (particularly necessary in 
the U n band). The da ta were r educed using procedures 
described in detail bv ISteidel et alJ (I2003T) . 

Additional images in the U n band covering the cen¬ 
tral 5' by 7' region of the WHT field were obtained 
in 2001 May using the blue side of the Low Resolu- 
tion Imaging Spectrometer on the Keck I 10m telescope 
llMcCarthv et al.lll998t ISteidel et al.ll2004l) . These were 
scaled, variance weighted and combined with the WHT 
data to form the final U n image. 

Properties of the reduced and stacked images are sum¬ 
marized in Tabled 

2.2. Near-IR Imaging 

We observed the HS1700+64 field during two observ¬ 
ing runs, 2003 May and 2003 October, using the Palo- 
mar 5.1m Hale telescope and the Wide Field Infrared 
Camera fWIRC: IWilson et alJl20 03L The images were 
taken in the K-short ( K s ) filter, obtained in a sequence 
of 30-second integrations, moving the telescope in a pre¬ 
scribed dither pattern after each sequence of 4 expo¬ 
sures. The conditions during both observing runs were 
generally good, with seeing ranging between 0"55 and 
0799, but the ambient temperature was warm result¬ 
ing in relatively high K s band sky backgrounds (~ 12.4 
mag arcsec -2 ). The WIRC camera employs a Rockwell 
HgCdTe “Hawaii-2” 2k x 2k array, providing a field of 
view of 8'.5 x 8i5 with a spatial sampling of 0.249 arcsec 
per pixel. 

The WIRC data from each (typically 54 minute) dither 
sequence were reduced, registered, and stacked using a 
set of IDL scripts written by K. Bundy (private commu¬ 
nication). The resulting images were then registered and 
combined using IRAF tasks. The final K s image repre¬ 
sents a total of 11.03 hours of integration in the parts 
of the image that received the full integration time. The 
final mosaic was trimmed to include only the parts of the 
image receiving 50% or more of the total integration time 
(~ 90% of the trimmed image received the full exposure 
time). The image w as calibrated w ith reference to the 
near-IR standards of lPersson et al.1 ( 1998 1. 

2.3. Spitzer/IRAC Imaging 

Observations with the Infrared Array Camera (IRAC) 
were obtained in 2003 October during the “In-Orbit 
Checkout” (IOC) of Spitzer. The observin g strategy and 
data r eductions are described in detail bv lBarmbv et al.l 
<200l . The IRAC observations were designed to over¬ 
lap as much as possible with the ground-based K s and 
optical images, but some compromises were adopted to 
avoid known foreground galaxy clusters in the field, so 
that the IRAC pointings were biased toward the SW rel¬ 
ative to the ground-based data, which are centered on 
the bright QSO. As detailed below, the IRAC images in 
the 4 IRAC bands (see Table 0 completely overlap over 
only a relatively small fraction of the field observed in 
the K s band described above, with the smallest overlap 
in the IRAC channel 1 (3.6/im) and channel 3 (5.8 /im) 
images, but much larger overlap in channels 2 (4.5/im) 
and 4 (8.0/rm). Figure ^ illustrates the regions covered 


in the IRAC images compared to those covered in the 
Palomar K s image. 

The photometric calibrations for the IRA C images 
were as determined by the IRAC team (iBarmbv et alJ 
l200l : however, we used a significantly different method 
for performing the photometry and estimating the pho¬ 
tometric uncertainty (see 13.21) . resulting in significantly 
more conservative (i.e. larger) estimated errors at a given 
flux density. 

3. PHOTOMETRY 
3.1. Optical and Near-IR 

The observables from our ground-based images are 72 
magnitudes and isophotal 7 Z — K s , G — IZ , and U n — G 
colors, with isophotes defined in the IZ detection cat¬ 
alog. In order to mimic the actual observing process, 
we therefore estimated uncertainties in isophotal color 
rather than magnitude for all optical/IR bands except 
IZ. The uncertainties were estimated separately for the 
optical and IR photometry, by adding large numbers of 
simulated galaxies of known magnitudes to the images 
and comparing their recovered photometry with their 
input parameters. The simulated galaxies were drawn 
from templates of exponential disks with scale lengths 
ranging from 0705 to 074, smoothed to match the seeing 
of the individual images. For the IZ — K s uncertain¬ 
ties, we selected galaxies at random within the range 
21.0 <7 Z < 26.0 and 1.0 <72. — I\ s < 5.5 and placed 
them at random positions on the IZ and K s images, 200 
at a time (to avoid overcrowding) for a total of 250,000 
fake galaxies. We then performed photometry as usual 
on the images, and compared the catalog of detections 
with the input list of simulated galaxies. We considered 
the quantity A[72.- K s \ = (IZ- K s ) meas (IZ A s)truei 
and binned in steps of 0.5 mag in 72 mea s and 0.2 mag 
in (IZ — K s ) me as- The mean value of A [IZ — K s ] is an 
estimate of the bias in the photometry; for typical ob¬ 
jects, A[72 — I\ s ] < 0.05 mag. The standard deviation 
er(A[72 — K s ]) provides an estimate of the IZ — K s uncer¬ 
tainty for objects that fall in each bin; this was typically 
er(A[72 — K s \) ~ 0.25 mag. IZ — K S uncertainties for each 
object are given in Tabled Uncertainties in 72, G — IZ 
and U n — G were determined in a similar fashion, by 
adding large numbers of fake galaxies to the U n , G and 
IZ images and comparing their recovered and input prop¬ 
erties, with the additional restriction that only those ob¬ 
jects whose recovered photometry met our selection crite¬ 
ria were considered. We binned in 72 and G — 72 to deter¬ 
mine the 72 and G—IZ uncertainties, which were typically 
a( A72) ~ 0.14 mag and <j(A[G —72]) ~ 0.07 mag, and in 
G and U n — G to determine U n — G uncertainties, typi¬ 
cally a(A[U n — G\) ~ 0.13 mag. Again uncertainties for 
each object are given in Table [2 Similar techniques for 
estimatin g photometric un c ertain ti es have been u sed in , 
e.g.. lAde lberger fc Steidell (120001) : ISteidel et alJ ll 2 00.31) : 
IShanlev et al .1 (120011) . 

3.2. IRAC 

We measured IRAC fluxes by first computing an em¬ 
pirical point spread function (PSF) for each channel by 
median averaging the signal from several isolated point 
sources located throughout the mosaiced IRAC data in 
the HS1700 field and normalizing the resulting PSF to 
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Fig. 1.— A schematic map of the survey field considered. The grayscale image is the Palomar 5m WIRC K s image, with the positions 
of objects with 1.4 < z < 2.8 marked with black circles. The colored contours (blue=IRAC/Ch 1; green=IRAC/Ch 2; yellow=IRAC/Ch 3; 
red=IRAC/Ch 4) indicate regions of the image that received > 10% of the total IRAC integrations time in each IRAC band, i.e., > 1.1 
hours. For full resolution figure, go to http://astron.berkeley.edu/~ aes 


have unit flux. We obtained accurate object positions 
( ;$ 0.2" rms uncertainty) for >95% of IRAC sources in 
the field from their AVband counterparts in the WIRC 
imaging of this held, or from their "A-band counterparts 
for the small fraction of objects detected in 1Z but not 
in K s . We then used the empirically-determined IRAC 
PSF to simultaneously model the emission of all objects 
within a sub-image of size 24" surrounding each object of 
interest, similar to methods generally used for photome¬ 
try of stellar objects in crowded fields. This procedure to 
some extent mitigates the effects of confusion noise, since 
the object positions are known a priori from the much 
higher spatial resolution K s band image, and allows for 
some de-blending of objects that are only partially re¬ 
solved in the IRAC i mages. Th e method is also similar 
to th at emp loyed by iFernandez-Soto et al.1 (1 9991 1 and 
IPapovich et alJ (1200111 for comparing ground-based near- 


IR and HST/WFPC-2 photometry. Photometric errors 
were computed by adding in quadrature the error in the 
source flux (Poisson noise) and the dispersion in mea¬ 
sured flux values found by fitting PSFs to 100 random 
positions containing no obvious sources near the object 
of interest. This latter estimate should account to a large 
extent for noise due to source confusion, which is partic¬ 
ularly relevant in the 3.6 and 4.5 /im bands. Except in 
cases of very bright sources, the background dispersion 
dominates the errors. The PSF full-width at half maxi¬ 
mum (FWHM) and approximate 5cr sensitivity limits for 
each channel are presented in Table 0 We omitted any 
measured fluxes with uncertainties larger than 0.5 magni¬ 
tudes from consideration in fitting the SEDs as described 
in s im In addition, we set the errors to 0.10 magnitudes 
for measurements whose formal uncertainties calculated 
using the above prescription were < 0.10 magnitudes, in 
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Redshift 


Fig. 2.— Redshift histogram for 72 spectroscopically iden¬ 
tified galaxies in the HS1700+643 field with near-IR and IRAC 
detections. This histogr am contains a signi ficant spike at z = 
2.30 ± 0.015, discussed bv ISteidel et alJ J2005D . 


order to account for systematic uncertainties in the pho¬ 
tometry that are not otherwise accounted for (e.g., source 
confusion on scales smaller than the PSF, estimation of 
local background, etc.). 

3.3. Sample Selection 

Because the redshift of HS1700+643 is z em = 2.72 
and we were interested in the properties of the galaxies 
in the “foreground,” the spectroscopic observations tar¬ 
geted photometrically selected galaxies that would neatly 
cover the redshift range 1.9 ;$ z ^ 2.7. To accomplish 
this we selected a combination of obj ects satisfy in g ei¬ 
ther th e “BX” photometric criteria of lAdelberger et~aTI 
(l2004bfl (z = 2.20 ± 0.32) o r the “MD” photometric cri¬ 
teria of lSteidel et alJ <1200311 (z — 2.79 ± 0.27). A total of 
100 objects (81 “BX” and 19 “MD”) has measured spec¬ 
troscopic redshifts z > 1.4 from observations with the 
Low Resol utio n Imaging S pectrometer, described in de¬ 
tail in lSteidel et alJ <1200411 . Because the surface density 
of BX candidates is ~ 6 times higher than that of MD 
candidates to the same 1Z = 25.5 apparent magnitude 
limit, the expected redshift distribution for a BX/MD 
photometric sample, based on the total of 1133 spec¬ 
troscopic redshifts for BX/MD candidates in all survey 
fields, is (z) = 2.24 ± 0.36. 

The sample of galaxies discussed in the rest of the 
paper is culled from the region of overlap between the 
ground-based U n GTZ images, the ground-based K s im¬ 
age, the 4 IRAC bands (see figure [TJ and the spectro¬ 
scopic observations of UV-selected galaxies. Within the 
8(3 by 7(8 region covered by the Palomar I\ s image, 
there is a total of 508 BX and MD-type galaxy candi¬ 
dates identified from the o ptical photometry using th e se- 
lection criter i a defin ed bv lAdelberger et a,U ll2004hH and 
iSteidel et alJ <1200311 . respectively; of those, 389 (77%) of 
the objects are significantly detected in the A' s -band im¬ 
age. Within the 64.7 square arcmin region covered by the 
K s image, we obtained spectra for 162 of the 508 pho¬ 
tometric candidates and spectroscopically identified 98 


objects. The identified fraction represents a 60.5% spec- 
troscop ic success rate, sim ilar to that of our z ~ 2 survey 
overal l (ISteidel et al.l2004H . As discussed in lSteidel et all 
<2001 . we expect the unidentifed objects have the same 
redshift distribution as the spectroscopic successes. Of 
the 98 objects with spectroscopic identifications, 92 are 
detected in the I\ s image to K s (Ve ga) ~ 22; of the 92 ob¬ 
jects with spectroscopic IDs and K s detections, 79 (86%) 
are at z > 1.4 and have a mean redshift and dispersion 
of 2.26 ± 0.30, entirely consistent with the expected red¬ 
shift distribution for all BX/MD candidates as discussed 
above. Of the 79 spectroscopically identified z ~ 2 galax¬ 
ies, 72 are detected in one or more of the IR AC passbands 
and included in the stellar population analysis. Invari¬ 
ably, objects detected significantly at 2.15 pn i are also 
detected in the IRAC 3.6 and 4.5 pm bands; the 7 galax¬ 
ies detected at 2.15pm but not with IRAC simply fall 
outside of the well-exposed regions of the IRAC point¬ 
ings. 

In summary, the 72 2 ~ 2 galaxies discussed in this 
paper are representative of a UV-selected sample to 
1Z = 25.5, but are objects which are also detected in 
both the near and rnid-IR. The only possible bias of this 
sample relative to any sample selected using the same 
UV criteria is that it excludes the ~ 23% of BX/MD 
candidates in the field that are not detected to K s ~ 22. 
The redshift distribution for this spectroscopic sample 
with near and mid-IR detections is shown in figure [2 
Very strikingly, there is a significant galaxy overdensity 
in this distribution at z = 2.30, which allows for an exam¬ 
ination of gal axy properties as a function of large-scale 
environment (ISteidel et alJl2005 L 

It is of interest, for some purposes, to ask what fraction 
of a mid-IR selected galaxy sample at the same redshift 
is represented by objects selected using the UV color cri¬ 
teria. This is a difficult question to answer, as there is 
as yet no large spectroscopic sample of mid-IR selected 
high redshift galaxies (and, as we argue below, photo¬ 
metric redshifts may be an inadequate approach to the 
problem). However, we can make a ro ugh estimate as 
follows. Recently, iBarmbv et al.l (2004) suggested a re¬ 
gion in IRAC 3.6/4.5/5.8 pm color space that may make 
it possible to select high redshift galaxies from IRAC col¬ 
ors alone, based on models and the observed locations of 
z ~ 3 UV-selected galaxies in IRAC color space. In fig- 
ure|3]we show the locations in the U n G7Z color plane of all 
of the objects detected in the I\ s band that are within the 
region of the IRAC images covered by IRAC at 3.6, 4.5, 
and 5.8pm and which have IRAC uncertainties of < 0.3 
mag in each band. 4 . Out of a total of 618 galaxies with 
K s , 3.6/im, 4.5/rm, and 5.8pm measurements, 443 are in¬ 
cluded in the optical photometric catalog to TZ = 25.5, 9 9 
of whi ch satisfy the criteria proposed by IBarmbv et alJ 
( 2004 *1 in the 3.6-5.8/rm color space. As shown in fig¬ 
ure □ 46 of the 99 sources satsifying the IBarmbv et alJ 
color criteria do indeed lie in the regions of U n G1Z color 
space from which z ~ 2 BM (z = 1.70 ± 0.34), BX 
(z = 2.20 ± 0.32], MD (z = 2.72 ± 0.3) and C/D/M 
(z = 3.0 ± 0.3) (ISteidel et al.112003 *) galaxies are drawn. 


4 The size of this region is approximately 30 arcmin 2 , although 
only ~ 40% of the objects in a iCs-selected catalog within this 
region are included, primarily because of the limited depth of the 
5.8/xm image. 
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Fig. 3. — Plot showing the location in optical U n GlZ color space of all objects in a K-selected catalog to K ~ 22. The regions within 
the solid black lines denote those used to select the galaxies discussed in this paper (“BX” and “MD”); the dotted curves denote regions 
containing “BM” (z = 1.7 =L 0.3, lower left) and “C/D” LBGs (z = 3.0 d= 0.3, upper left). Objects satisfying the “high z” galaxy criteria 
suggested by Barmby et al. (2004) based on IRAC colors are marked in red. As discussed in the text, a significant fraction (~ 40%) of the 
IRAC “high z” objects in the optical catalog appear to lie at 2 : < 1.4. 


However, there is a significant number of galaxies (~ 40) 
which satisfy the Barmby et al. criteria but are located 
in regions inhabited primarily by z ~ 0.7 — 1.5 star form¬ 
ing or passive galaxies in optical color space (see, e.g., 
lAdelberger et al.ll2004bll . Another ~ 10 have optical col¬ 
ors that are outside of the UV selection windows, but are 
consistent with more reddened sta r forming galaxies at 
z ~ 1.5- 2.5 (cf. lDaddi et al.ll2004all . Of the 175 galax¬ 
ies that are well detected in the K s , 3.6, 4.5, and 5.8 
(im bands but a bsent from the optical catalog, 81 sat¬ 
isfy the iBarmhv et a,l.l color criteria. Assuming that the 
same percentage (~ 40%) of these are likely to be in the 
foreground as for the optically detected galaxies, ~ 50 
galaxies to K s ~ 22 are plausibly at z ^ 1.5 s . Most of 

5 29 of the 81 objects are actually detected in the 1Z band to 
TZ ~ 26 but are not included in figure [3] because the colors are not 


the galaxies in this red sample are relatively faint in the 
K s band, with only 12 of 81 having K s < 20, of which 
perhaps 6-8 are high redshift objects. 

In summary, of the 180 galaxies well detected in the 
in t he A,. 3.6, 4.5, and 5.8 (im bands that also satisfy 
the iBarmbv et al.l criteria, 46 would be selected based 
on U n GlZ colors to be at z £ 1.5, ~ 70 (~ 40%) are 
likely to be at z ^ 1.5, while the remaining ~ 60 galaxies 
are plausible z £ 1.5 candidates that do not satisfy the 
U n GlZ selection criteria. Thus, we conclude that ~ 50% 
of galaxies at z ~ 2.3 ± 0.4 to an apparent magnitude 
limit of K ~ 22 (K ~ 23.8 on the AB system) would 
be selected using the rest-UV color criteria to 7 Z = 25.5. 
This statistic is broadly consisten t with the results of the 
FIRES survey llFranx 7Tanil2003h which targeted galax- 

reliably measured. 
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ies with colors J — K > 2.3, most of which are similarly 
optically faint, as well as with the “IRAC EROs” dis¬ 
cussed bv lYan et al.l ( 200 1 i. We defer further discussion 
of the overlap and distinctiveness of samples selected in 
various ways to future work. 

4. POPULATION SYNTHESIS MODELING 
4.1. Sample 

In order to explore the stellar populations of UV- 
selected z ~ 2 galaxies in the HS1700+64 IRAC pointing, 
we considered only the subsample of 72 galaxies with 
measured redsliifts, K s magnitudes, and IRAC detec¬ 
tions in at least one channel. As discussed below, the 
set of colors and magnitudes for BX and MD galaxies is 
not sufficient to simultaneously determine accurate pho¬ 
tometric redshifts and stellar masses for galaxies that do 
not have spectroscopic redshifts. 

The photometric measurements used to model the stel¬ 
lar populations are given in Table 2. Note that all of the 
magnitudes and colors are on the AB system, with the 
exception of the 1Z — K s color, for which the K s mag is 
on the Vega system for easier comparison to other work 
on similar high redshift galaxies. 6 For the 13 galaxies at 
z > 2.5, we did not include the U n — G color in the SED 
fit, since absorption by H I in the intergalactic medium 
begins to have a strong effect on the measured XJ n mag¬ 
nitude at these redshifts. 

Before modeling the galaxy magnitudes and colors, we 
corrected the U n or G magnitude for the effects of galac¬ 
tic Lya emission and absorption. Lya equivalent widths 
were measured from the LRIS spectra used to identify 
redshifts. These corrections were applied to galaxies at 
z = 1.7 —1.9 (U n ) or z > 2.47 (G). We also have Ha and 
[Nil] emission line measurements for four galaxies in the 
HS170 0 sample (BX691, BX717, MD103, and MD109, all 
from lErb et al.1 l|2003l) l. We corrected the observed K s 
by 0.10, 0.17, 0.08, and 0.10 magnitudes, respectively to 
account for the contribution of nebular emission to the 
observed flux. K s magnitudes for the remainder of the 
galaxies in the redshift range 2.016 < z < 2.504 remain 
uncorrected. Based on a much larger sample of Ha spec¬ 
tra of galaxies in the same redshift ran ge (which will be 
discussed in detail bv lErb et al. 1, 2005 )). the typical cor¬ 
rection to the broad-band K s flux is +0.09 ± 0.07 magni¬ 
tudes, although some of the most extreme cases have Ha 
equivalent widths large enough to change the K s magni¬ 
tude by as much as several tenths of a mag. The possible 
impact of the K s corrections on the modeling results is 
discussed below. 

Only 22% of the BX/MD photometric candidates with 
both K s photometry and measurements in at least one 
IRAC band are spectroscopically identified, and we re¬ 
strict the stellar population modeling to this spectro¬ 
scopic sample. To test the potential success of model¬ 
ing photometric candidates without redshifts, we tried 
to fit simultaneously the redshift and stellar popula¬ 
tion parameters for the sample of galaxies with spec¬ 
troscopic redshifts, pretending that the redshift was 
an undetermined free parameter. With only U n GlZK s 
plus IRAC photometry, we obtained an accuracy of 
(| ^ph,ot z spec 1/(1 + z spec )) — 0.09. At the mean red- 

6 K s magnitudes or 7 Z — K s colors may be converted to the AB 
system by adding 1.82 to the tabulated values. 


shift of the sample, this error corresponds to a typical 
|Zphot — z spec | of 0.3, while in several cases the redshift 
estimates fail catastrophically, with \z p h. 0 t — z spec \ > 1. 
We note that this accuracy is very similar to that found 
for photometric redshift estimates in the FIRES sur¬ 
vey, based on 39 galaxies with spectros copic redshifts 
in the Hubble Deep Field South (HDF-S) (Rudnic k et alJ 
12003 ; lLabbe et al.1 ' 2003 1. Perhaps more significantly, 
the stellar population parameters inferred when simul¬ 
taneously fitting for the redshift do not uniformly agree 
with those found when the redshift is known. Specifi¬ 
cally, we find a typical fractional stellar mass difference 
of (\M* hot — M* pec |/M* pec ) = 0.56 + 1.23, which is larger 
than any of the fractional stellar mass uncertainties we 
find due to photometric uncertainties, or uncertainties in 
our choice of star-formation history for single-component 
models (see discussion in section E3 - Not only is the 
scatter in this fractional offset quite large, but the distri¬ 
bution is skewed toward overestimating the mass, such 
that the average ((M* hot — M* pec )/M* pec ) = 0.37+1.30. 7 
This exercise demonstrates that, with only U n GTZK s plus 
IRAC photometry, even for objects known to have red¬ 
shifts in the targeted range (which presumably reduces 
the probability of catastrophic failures), considerable un¬ 
certainties would be introduced if objects without spec¬ 
troscopic redshifts were included. 

4.2. Modeling Procedure 

Foll owing a pr ocedure simila r t o that of lShanlev et alJ 
Ml and iShanlev et all IMMJ). we apply 
iBrnziia.l V Chariot! ( 2003 1 models to fit the SEDs 
of the UV-selected z ~ 2 galaxies in the HS1700 
Spitzer pointing. We used models with solar metallicity 
and a Salpeter initial mass function (IMF) extend- 
ing from 0.1 — lO O Mm. Also, as recommended by 
fBruzual fc Chariot! ( 2003 *1. we ad opted the Padova 
1994 stellar evolution tracks. iPa.novich et all ( 12001 h 
have investigated how the best-fit parameters depend 
systematically on the choice of metallicity and IMF, 
though we will not include such a discussion here. We do 
note that, at least for galaxies bright in the rest-frame 
optical (K s < 20), solar metallicity appears to be a 
good approximation, and typical z ~ 2 galaxies appear 
to ha ve meta llicitie s that are only slightly lower than 
solar t Shaolev et al . 200 41. We also note that it is well 
established that a Salpeter IMF over-predicts the num¬ 
ber of stars less massive than 1 M 0 (or, alternatively, 
that it over-predicts the stellar M / L) compared t o 
observations in the local universe (e.g., [ Bell et al . I®. 
We use the Salp eter IMF to facil i tate comparison wit h 
other wor k (e.g..lCole et al .1120011 iShaolev et al .1 l200Hl : 
use of the iChahried ^^^TlMF, the IMF proposed by 
iBaldrv fc Glazebrook 120031 , or the “diet” Salpeter IMF 
used bv iBell et all (l2003^h would reduce the inferred 
stellar masses by a factor of ~ 1.5 — 1.8. D ust ex ti nction 
is taken into account with a iCalzetti et al.l ( 1200011 
starburst attenuation law. iR.eddv fc Steidell (I2004H have 
shown that using the Calzetti law and continuous star- 
formation models to infer unobscured star-formation 

7 When we exclude five galaxies with wildly discrepant M% hot 
relative to M* pec (all are overestimates), we find a relatively sym¬ 
metric distribution around zero, with ((M* hot — M* pec )/M* pec ) = 
0.07 + 0.44. 

















































rates from UV colors and magnitudes accurately pre¬ 
dicts the average X-ray and radio continuum fluxes of 
2 ~ 2 star-forming galaxies, though the object to object 
scatter may be large. In order to determine how well 
the stellar population parameters can be constrained, 
we investigated a range of star-formation histories of 
the form SFR(t) oc exp(—i/r), with e-folding times 
of r = 0.01, 0.05, 0.10, 0.2, 0.5, 1, 2, and 5 Gyr, as 
well as continuous star-formation (CSF) models. For 
the main analysis, we did not consider more complex 
star-formation histories, such as short-duration bursts 
occurring at random past times superimposed on a 
smoothly declining star-formation history. The data for 
most galaxies do not place strong constraints among 
the range of simple, single-component star-formation 
histories considered, and therefore the use of a more 
complex model does not seem justified. We only br iefly 
employ extreme two-component models in section IOI 
in order to determine how much stellar mass could be 
contained in these galaxies in maximally old underlying 
bursts. 

Our goal was to constrain the stellar population pa¬ 
rameters of each galaxy. The fitted parameters included: 
dust extinction (parameterized by E(B — V)), age (t s f), 
star-formation rate (SFR), stellar mass (M*), and star- 
formation history ( t). F or each r, we used a grid of 
IBruzual fc Chariot (2003) models with ages ranging be¬ 
tween 1 Myr and the age of the Universe at the redshift 
of the galaxy being modeled and extinctions ranging be¬ 
tween E(B — V) = 0.0 and E(B — V) = 0.7 (though none 
of the galaxies has a best-fit E(B — V) >0.5). For each 
combination of E(B — V), t B f , and r, we computed the 
full U n GlZK s plus IRAC 3.6 — 8.0/iin SED of the model 
placed at the redshift of the observed galaxy. In com¬ 
puting this SED, we reddened the intrinsic model galaxy 
spectrum by dust, and further attenuated it in a man¬ 
ner simulating abs orption by the intergalactic medium of 
neutral hydrogen iMadaulll995th which affects the pre- 
dicted U n and G magnitudes. IBruzual fc Chariot ( 1200 3) 
CSF models are by default normalized to have star- 
formation rates of lM 0 yr — , while declining models are 
normalized to have a total mass of lMg at t = oo. We 
searched for the model normalization that minimized the 
value of x 2 with respect to the measured SED, which 
yielded an estimate of the star-formation rate and stellar 
mass. The combination of E(B — V), t s f, r and nor¬ 
malization yielding the lowest y 2 was called the overall 
“best-fit”. Figure 0 shows the overall best-fit model to¬ 
gether with the observed data points for each of the 72 
galaxies in the sample. In most cases the “best-fit” r 
model is not significantly better than that of the best¬ 
fitting model for every other value of r considered, as 
discussed below. 

4.3. Systematic Uncertainties 

Even neglecting the effects of photometric uncertain¬ 
ties, parameters derived from SED fitting are sub¬ 
ject to significant degeneracies and systematics, pri¬ 
marily because of the inability to constrain the star- 
forma tion history I Sh aolev et all 120011: iPanovich et al.l 
12001 . These uncertainties apply even using the simplest 
single-component models, as we do here, since the best-fit 
age, extinction, and star-formation rate all depend sys¬ 
tematically on the t used to describe the star-formation 


history. Stellar mass estimates are more robust in the 
face of uncertainties in the star-formation history. Here 
we review the systematic uncertainties in age, extinc¬ 
tion, star-formation rate, and stellar mass that result 
from population synthesis modeling of galaxy SEDs at 
high redshift. 

The inferred age, t s f, for a given value of r is con¬ 
strained mainly by the magnitude of the age-sensitive 
Balrner break. A given Balmer break strength indicates 
the relative number of A and O stars and corresponds 
to older ages when larger values are assumed for t (i.e. 
a given relative number of A/O stars is reached at later 
times for more gradually declining star-formation histo¬ 
ries). 

The inferred extinction, E(B — V), also depends sys¬ 
tematically on the star-formation history: declining star- 
formation histories with t s f /r > 1 with less dust extinc¬ 
tion can produce similar G—1Z and 1Z—K S colors to those 
of a model with a larger r, smaller t s f /r and more dust 
extinction. In this case, the reddening of the rest-frame 
UV slope is caused by a mixture of stars with later spec¬ 
tral types on average for models with smaller r, whereas 
dust extinction accounts for the UV reddening in contin¬ 
uous star-formation models or those with larger values 
of r. 

The systematic dependence of the inferred SFR on the 
value of r that parameterizes the star-formation history 
stems from two effects. First, the star-formation rate 
is derived from the extinction-corrected UV-luminosity. 
It therefore depends on the best-fit value of E(B — V ), 
which as we discussed above, depends on the star- 
formation history used to model the observed SED. 
Second, a conversion from extinction-corrected UV- 
luminosity to SFR is required. This conversion depends 
on both t s f and r, and for a given galaxy, is gener¬ 
ally an increasing function of r. The conversion from 
UV-luminosity to SFR equilibr ates after a b out 100 Myr 
for models with r > 1 Gyr llKennicu How¬ 

ever, the ratio of SFR to UV luminosity declines as 
a function of time for models with smaller r. In the 
range of parameter space spanned by the galaxies in 
our sample, r = 100 Myr models typically have UV-to- 
SFR conversions 1.2 — 1.8 times lower than continuous 
SFR models that fit the same colors. While not likely, 
based on external (spectroscopic) information, models 
with r = 10 Myr have much more extreme conversions 
from UV-luminosity to SFR, as much as 25 — 1000 times 
lower than the limit approached by the models with large 
r. This is easily understood since the largest r = 10 Myr 
age that can reasonably fit the observed colors of the 
galaxies in our sample is roughly t s f ~ 100 Myr. At 
this point, the star-formation rate is roughly 20000 times 
lower than at t s f = 0. However, the UV-luminosity has 
only declined by a factor of ~ 20, because while there 
are no O-stars left, the lifetime of late B-stars is roughly 
100 Myr, so they have not yet disappeared. 

Our single-component-model estimates of stellar 
masses are not plagued by the same amount of system¬ 
atic uncertainty, especially when considering more real¬ 
istic exponential timescales, r > 100 Myr. There has 
been much discussion in the literature about the value 
of rest-f rame near-IR photomet r y in measuring stellar 
masses llClazebrook et al.1120041 iBell et alJ 1200311 . In 
fact, as shown by the model fits to BX490 and BX505, 
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Fig. 4. — Examples of confidence intervals obtained from Monte Carlo modeling for three different galaxies selected to span a range of 
inferred properties, as described in the text (BX879, BX691. and BX563). The confidence intervals have been projected onto the spaces 
of (E(B-V),Age) and (Age,Mass), and the dark and light contours correspond to 68% and 95% confidence ranges, respectively. While the 
time constant, r, is treated as a free parameter ranging from 10 Myr to oo, the grid in r contains only ten distinct values. The discrete 
nature of r translates into discrete regions of parameter space, that would be connected smoothly if r were allowed to vary continuously. 
This effect is conveyed most clearly in the (E(B-V),Age) plots, since both parameters systematically depend on r. For full resolution figure, 
go to http://astron.berkeley.edu/~aes 


which have very similar rest-frame near-IR magnitudes, 
rest-frame near-IR photometry alone is insufficient to 
constrain the stellar mass-to-light ratio of star-forming 
galaxies to better than a factor of 10 (see also M5.4H . 
However, with rest-frame UV, optical, and near-IR pho¬ 
tometry, the near-IR M/L ratio (and therefore the stel¬ 
lar mass) is reasonably well-determined, almost indepen¬ 
dent of the value of r used to fit the SED. With mea¬ 
surements of the rest-frame UV slope (observed G — 1Z 
color at high redshift), the rest-frame UV/optical color 
(observed 1Z — I \ s ), and rest-frame UV/near-IR color 
(observed 1Z — IRAC), the variance in inferred near-IR 
M/L ratio due to uncertainties in the form of the star- 
formation history (for r > 100 Myr) is typically 10%. 
The relationship between (G — TZ/lZ — K s , 7 Z — IRAC) 
and M/L ratio starts to break down for star-formation 
histories with r < 100 Myr. The systematic uncertainty 
in the AL/L ratio therefore increases, except in the cases 
where t s f < 10 Myr (and < r for all star-formation histo¬ 


ries under consideration), when all star-formation histo¬ 
ries yield the same best-fit stellar population parameters. 
However, for most galaxies in the sample, there is no indi¬ 
cation that models with r < 100 Myr with t s f /r > 1 are 
preferred, and in some cases such star-formation histories 
are actually ruled out (see discussion in section l(7~21l . 

External constraints can be used to evaluate the likeli¬ 
hood that these extreme “decaying burs t” sta r-for mation 
histor ies are realistic. For example, iReddv fc Steidell 
ll200l have shown that the mean extinction-corrected 
star-formation rate for z ~ 2 galaxies in the Chandra 
Deep Field North, inferred from rest-frame UV luminosi¬ 
ties and colors, assuming > 100 Myr of continuous star 
formation and the Calzetti extinction law, agrees very 
well with the average star-formation rate derived from 
stacked X-ray and radio images of the same galaxies. If 
r < 100 had been assumed, the average UV-derived star- 
formation rate would have been several times smaller, 
and no longer in agreement with the star-formation rate 
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inferred from the average X-ray or radio continuum flux. 
In principle, the degeneracies that affect primarily the in¬ 
terpretation of the UV continuum slope can be removed 
if there is external information on the number of massive 
stars still present in the galaxies. For example, using 
observed Ha fluxes and/or high quality far-UV spectra, 
it is relatively easy to establish whether O stars are still 
present in the galaxy spectra. In general, the high success 
rate in detecting Ha emission from UV-selected galax¬ 
ies, and in particular the fact that the star-formation 
rate inferred from Ha fluxes is never smaller th an th at 
inferred from the UV continuum fe.g.. lErb et a.l,112003 1. 
suggests that it must be quite rare to catch a starburst 
with tsF/T » 1 and r < 100 Myr. Additionally, in 
the individual LRIS rest-frame UV spectra with suffi¬ 
cient S/N, one sees evidence for the existence of massive 
stars from the presence of C IV and Si IV P-Cygni stel¬ 
lar wind features, and, in some cases Lya emission. If 
t < 100 Myr and t s f /t were much greater than 1 in most 
cases, we would not expect to find such strong evidence 
for massive O stars in the rest-frame UV spectra of these 
objects. 

The systematic effects discussed above do not include 
the uncertainties that result from including models with 
underlying short bursts of 
at random times in the p 
We did not add that layer of complexity to our mod¬ 
eling procedure, given the difficulty we found in dis¬ 
tinguishing among even simple exponentially declining 
and continuous star-f ormation histories. However, as 
iPanovich et al.l ( 200 11 demonstrated, underlying maxi¬ 
mally old bursts with t s f » r can increase the M/L by 
a factor of a several without having an appr eciable e ffect 
on the rest-frame UV-to-optical colors. Clay.e hrook et all 
lt20f)l find that including random bursts into the range 
of star-formation histories used to fit the VIz'K colors 
of 1.0 < ^ < 2.0 galaxies causes the derived masses to in¬ 
crease by less than a factor of two, although the amount 
by which one might be underestimating the true stellar 
mass will depend on the observed galaxy colors. We re¬ 
turn to this question in section IQ 


diff erent strengths oc cu rring 


4.4. Monte Carlo Confidence Intervals based on 
Photometric Errors 

In addition to the best-fit stellar population parame¬ 
ters for each galaxy, we also found the parameter confi¬ 
dence intervals allowed by the photometric uncertainties. 
A large set (n = 1000) of fake SEDs was generated for 
each galaxy, by perturbing the observed SED in a man¬ 
ner consistent with the estimated photometric error bars 
(sections Id.II and Id.211 . Each perturbed SED was fit in 
exactly the same manner as the observed SED, and a set 
of best-fit stellar population parameters was chosen to 
minimize y 2 with respect to the perturbed SED. 

We determined the confidence intervals over E(B — V ), 
t s f, SFR, M*, and r, such that 68.3% and 95.4% con¬ 
fidence intervals contain the corresponding fractions of 
the Monte Carlo realizations. More than half of the 
galaxies in our sample have 95.4% confidence intervals 
that span all r values, indicating the difficulty of dis¬ 
tinguishing among the different star-formation histories 
with broad-band photometry alone. In practice, the most 
useful number we derive from the Monte Carlo confidence 
intervals for each galaxy is the uncertainty on the stellar 
mass. For this quantity, it is important to know not only 
the random effects of photometric errors, but also the 
systematic effect of allowing r to vary. The ler Monte 
Carlo stellar mass uncertainties in column 7 of Table 3, 
derived from simulations in which r was allowed to vary, 
reflect the uncertainty in r. We also list the Monte Carlo 
stellar mass uncertainties for CSF star-formation histo¬ 
ries in Table 3. 

To demonstrate the covariance of different parame¬ 
ters, and at least some of the range of stellar popula¬ 
tions present in our sample, we show in figure @1 the 
t s f vs. E{B — V) and M* vs. t s f confidence inter¬ 
vals derived from Monte Carlo simulations. These con¬ 
fidence intervals do not assume a fixed value for r, but 
are the results of simulations in which r was allowed to 
vary. The confidence intervals displayed include those of 
BX879 (z = 2.308), which appears to be a very young 
( t s f < 10 Myr), dusty stellar population with high in¬ 
stantaneous star-formation rate and low stellar mass, all 
of which parameters are independent of r because of the 
extremely young age; BX691 (z = 2.190), an evolved 
stellar population with t s f/r > 1, r constrained to be 
between 0.5 and 2 Gyr, moderate current star-formation 
rate and large stellar mass; and BX563 (z = 2.296), for 
which there are no constraints on the star-formation his¬ 
tory, but the stellar mass of which is close to the sample 
median. The large extent of the BX563 confidence inter¬ 
val in t s f vs. E(B — V ) space demonstrates some of the 
fundamental limitations of population synthesis model¬ 
ing. 


4.5. Caveats to the Model Results 

The goal of population synthesis modeling is to in¬ 
fer physical quantities from observed galaxy magnitudes 
and colors, in particular the stellar mass. In the best- 
case scenario we would also like to place constraints 
on the timescale of the current star-formation episode, 
and therefore the level of dust extinction and the dust- 
corrected star-formation rate. There are certain limita¬ 
tions of the modeling procedure, however, that are worth 
pointing out, in addition to the systematic uncertainties 
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Fig. 5.— Observed and model SEDs for all 72 of the galaxies with spectroscopic redshifts and measurements in the K s band and at 
least one IRAC channel, sorted by redshift from lowest to highest. The numbers for the inferred stellar mass are taken from the best-fitting 
exponentially declining model, which is shown shown with a solid curve. The dotted curve represents the best-fit model when the IRAC 
data points are excluded. 


discussed in section FOl 

First, it is not always possible to find a statistically ac¬ 
ceptable fit to the observed galaxy colors, especially using 
relatively simple models. Seven galaxies in the sample 
have best-fit models with significantly higher y 2 values 
than those of the rest of the sample. Two of these galax¬ 
ies (BX490 (z = 2.403) and BX794 (z = 2.253)) have 
large 1Z — K s residuals that are likely to be explained 
in terms of contamination from Ha emission in the I\ s 
band. Another galaxy, BX561 (z = 2.426), has a good 
fit except for the IRAC 8.0 /im point, which yields a sig¬ 
nificant positive residual, possibly due to the presence of 


an obscured AGN (see H3Ji. MD157 shows rest-frame 
UV spectroscopic evidence for hosting a broad-line AGN, 
which probably contaminates the K s flux (strong Ha) 
and IRAC 8.0 gm band. The SEDs of the remaining 
three galaxies (BX581 (z = 2.406), BX756 (z = 1.738), 
and BX681 (z = 1.740)) are very difficult to understand, 
with both large K s and IRAC residuals that are both 
positive and negative, and don’t stem from obvious ir¬ 
regularities in any of the imaging data. 

A second issue is the apparent discrepancy between 
model parameters based on the best fit to broad-band 
photometry and those inferred from rest-frame UV spec- 
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Fig. 5.— (continued) 


tra, which applies for a small number of galaxies in the 
sample: BX1032 (z = 1.883), BX1087 (z = 1.871), and 
BX536 (z = 1.977). Not only did we obtain statisti¬ 
cally acceptable fits to the observed photometry of these 
galaxies, but also the results of the Monte Carlo simu¬ 
lations appear to rule out star-formation histories with 
r > 200 Myr, and to favor best-fit t s f/r significantly 
greater than unity, with quite modest star-formation 
rates of a few M@ yr -1 . The broad-band photometry of 
these galaxies implies star-formation histories in which 
the bulk of stellar mass formed in the past, and the cur¬ 
rent star-formation rate is much lower than the past av¬ 
erage. We also have rest-frame UV Keck I/LRIS spectra 
for these galaxies, which can be used as a consistency 


check on the stellar population parameters derived from 
the SED-fitting procedure. The spectra of BX1032 and 
BX1087 both show prominent C IV P-Cygni features, 
which are produced in the winds of massive stars, associ¬ 
ated with active ongoing star formation. The rest-frame 
UV spectra of BX1087 and BX536 both contain broad 
He II emission, which is produced by Wolf-Rayet stars, 
another indication of a young stellar population. Finally, 
the spectrum of BX1032 exhibits nebular emission lines 
such as CIII] and He II, further indications of active cur¬ 
rent star formation. In short, the rest-frame UV spectra 
of these t s f/r >> 1 galaxies contain features that are 
not consistent with the best-fit models used to describe 
the rest-frame UV to near-IR SED, which is all the more 
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Fig. 5.— (continued) 


striking since the models appear fairly well-constrained. 

Clearly, there are cases in which the simple, single¬ 
component SED-fitting procedure fails to capture the 
true nature of the galaxy stellar population. Rest-frame 
UV and optical spectra, and X-ray, submillimeter, and 
radio fluxes, provide valuable complementary informa¬ 
tion that must be included to better-constrain all of 
the stellar population parameters. The limitations in¬ 
herent in SED fitting, without external constraints (and 
especially without spectroscopic identifications), should 
be borne in mind when considering the resulting stel¬ 
lar population parameters. On the other hand, the stel¬ 
lar masses inferred for these three galaxies are largely 
independent of the star-formation history assumed, and 


therefore we include them in all further discussion of stel¬ 
lar masses. 

5. RESULTS 

5.1. The Stellar populations of z ~ 2 Galaxies in the 
HS1700 Field 

In Figure E3 we show the observed optical, near-IR, 
and Spitzer photometry and best-fit models for the en¬ 
tire sample of UV-selected 2 ~ 2 galaxies in the HS1700 
field with spectroscopic redshifts, K s , and Spitzer/ IRAC 
photometry. The so l id lin es indicate best-fit exponential 
iBruzual & Chariot! l)2003l) models including the IRAC 
measurements, whereas the dotted lines indicate fits to 
only U n GRK s photometry. The stellar mass inferred 
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Fig. 5.— (continued) 


from the best-fit to the entire SED (including IRAC pho¬ 
tometry) is also shown, along with the redshift of each 
object. In Table 0 we list the best-fit stellar popula¬ 
tion parameters for each object, including the entire op¬ 
tical through IRAC SED. Parameters are listed for the 
best-fit model assuming r = oo (CSF), and also for the 
exponentially declining model that formally yielded the 
lowest value of y 2 when r was allowed to vary as a free 
parameter. There is a subset of galaxies with extremely 
young best-fit ages (t s f < 10 Myr), for which all val¬ 
ues of r yield nearly identical stellar population param¬ 
eters. These galaxies include: BX928, BX1075, BX794, 
AID 128, BX879, BX490, MD155, MD152, BX588, and 
AID 126. In such cases where t s f is smaller than any of 
the t values considered, the star-formation history (as 
parameterized by r) becomes irrelevant. 

From both Figure 0 and Table 0 it is clear that a 
large range of rest-frame UV to rest-frame near-IR SEDs 
is included in our sample. In particular, we highlight the 
variety of rest-frame UV/near-IR colors, which translates 
directly into a large range of inferred near-IR mass-to- 
light ratios. 

5.2. Distribution of Stellar Masses 

Figure El shows the distribution of stellar masses in¬ 
ferred from the optical-IRAC SEDs. The mass for each 
object is taken from column 7 of Table 0 The nearly 
log-normal distribution of stellar mass can be charac¬ 
terized by a mean and standard deviation logM* = 
10.32 ± 0.51 Mq. For the CSF models (the results of 
which are also compiled in table 0 , we find a distribu¬ 


tion logM^(CSF) = 10.30 ± 0.53. Figure 0 compares 
the stellar mass obtained assuming the best-fit r model 
to the stellar mass implied by the best-fit CSF models. 
The agreement, both in terms of the mass distribution 
and the masses of individual objects, is excellent, with 
no obvious systematic trend or offset. 

For some purposes, it is of interest to break down the 
inferred stellar mass distribution as a function of ob¬ 
served AT s -band magnitude, since a number of other sur¬ 
veys of galaxies at similar redshifts are selected in that 
band. The results are summarized in Table 0 where the 
K s limits are chosen to correspond to the limits imposed 
by some of these other surveys. As expected, there is 
a clear trend of increasing inferred stellar mass with in¬ 
creasing rest-frame optical luminosity, although, as dis¬ 
cussed below, there is significant scatter in rest-optical 
M/L at a fixed mass or luminosity. 

The typical (median) Alonte Carlo uncertainties in in¬ 
dividual stellar mass estimates are ~ 30 — 35%, derived 
from simulations in which r was allowed to vary as a free 
parameter. If we model all galaxies in the sample with 
CSF star-formation histories, we find that the median 
Alonte Carlo stellar mass uncertainty is ~ 35%. These 
values indicate that, for the same photometry and un¬ 
certainties, CSF models actually yield stellar mass un¬ 
certainties at least as large as the ones derived when r is 
allowed to vary as a free parameter. 

5.3. Results with and without Spitzer/IRAC Photometry 

One of the fundame ntal limitations of the ana lysis 
in lShaolev et all (I200H) and iPanovich et all 11200111 was 
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Fig. 6.— Histograms of best-fit stellar masses, with (shaded) 
and without (unshaded, dashed) IRAC photometry. As discussed 
in the text, the overall mass distribution and the total inferred 
stellar mass are similar to within ~ 15%, although the uncertainties 
in the measurement of individual stellar masses are smaller by a 
factor of ~ 1.5 — 2 when the IRAC data are included. 



Fig. 7. — A comparison of the inferred stellar masses for constant 
star formation (CSF) models and for the best-fit declining star for¬ 
mation models (best r) for the 72 galaxies with IRAC measure¬ 
ments. The error bars are as presented in Table HI and are based 
on Monte Carlo modeling given the photometric errors. Note that 
there are no obvious systematic differences between the CSF and r 
models, although in some cases the r models result in significantly 
better agreement with the observed data. 


the the inability to discriminate among different star- 
formation histories with only rest-frame UV and optical 
photometry. We have shown that, even with the addi¬ 
tion of Spitzer /IRAC photometry, this inability to con¬ 
strain star-formation histories persists for most galaxies 
in our z ~ 2 sample. The broad-band SEDs are par¬ 
ticularly degenerate with respect to age, dust content, 
and star-formation rates, while the inferred stellar mass 
suffers from much less systematic uncertainty. Another 
question raised in the pre-Spitzer era was whether or not 
rest-frame near-IR photometry of z > 2 galaxies would 
uncover considerably more stellar mass from old bursts 
than had been inferred from SEDs whose longest wave¬ 
length measurement was the ground-based K s band, cor¬ 
responding to the rest-frame optical at the redshifts of 
interest here. Such stellar composite populations might 
have been revealed through IRAC fluxes in excess of the 
rest-frame near-IR extension of model fits to the rest- 
frame UV through optical SEDs. 

Figure 0 shows the relationship between 
M*(No IRAC), the stellar mass inferred when the IRAC 
photometry is excluded from the fit, and M*(IRAC), 
the stellar mass derived using IRAC photometry. 
Figure El displays this relationship in one dimension, 
with the distribution of M*(No IRAC) indicated as an 
unshaded histogram in front of the shaded histogram 
of M*(IRAC). What both figures demonstrate is that, 
on average, the stellar mass inferred with and without 
IRAC data agree without a large systematic offset, 
though the scatter in this relationship is significant 
and larger than the typical (formal) uncertainties in 
the individual stellar mass estimates using the full 
SED. Table 0] compares these inferred stellar mass 
distributions as a function of limiting observed K s flux. 

The average of log M*(No IRAC) - log M*(IRAC) 
=0.06, with a standard deviation of 0.35 dex; the me¬ 
dian difference is only 0.01 dex. We find only five galax¬ 


ies whose optical-IRAC SEDs imply stellar masses more 
than a factor of two larger than that inferred from the 
optical-A'., SEDs. The most dramatic such example is 
BX898, whose observed 4.5/im flux is ~ 3.5 times larger 
than that predicted by the fit to the optical-A),- SED 
alone. The inferred stellar mass is ~ 12 times larger if 
the IRAC 4.5/im point is included in the SED fit. Un¬ 
fortunately, BX898 has a significant IRAC measurement 
only at 4.5/im, and contamination of this measurement 
(or, a statistical outlier in the observed K s measurement) 
cannot be ruled out. 

It is more common for the fit to the optical-A"), SED 
to overpredict the stellar mass relative to that inferred 
from the optical-IRAC SED. We believe that this dis¬ 
crepancy arises in most cases because of the contamina¬ 
tion of the K s continuum flux by Ha + [Nil] line emission 
at 2.0 < z < 2.5. The addition of the IRAC photome¬ 
try to the SEDs of galaxies with strong Ha emission 
sometimes reveals the K s data point as a positive outlier 
relative to the overall SED fit (e.g. BX794 and BX490), 
and the best-fit Balmer break is smaller than the ob¬ 
served 72. — K s color alone would imply. In the absence 
of IRAC photometry, the K s measurement strongly in¬ 
fluences the inferred stellar mass because it is the only 
proxy for the strength of the Balmer break, resulting in 
a bias toward somewhat higher values. In support of 
this interpretation, for the four galaxies in the HS1700 
A' S /IRAC pointings with published Ha measurements, 
we find an average ratio of 1.55 between M*(No IRAC) 
and M*(IRAC) when no correction for Ha emission is 
applied. When Ha emission is corrected for, we find an 
average M*(No IRAC)/M*(IRAC) of 1.11, i.e. closer to 
unity. More complete Ha surveys are necessary to re¬ 
move the K s emission line bias, especially in order to 
model objects for which I\ s is the reddest available band 











16 



Fig. 8. — A comparison of the inferred stellar masses for the 
best-fit declining star formation models (best r) for the 72 galaxies 
with IRAC measurements, derived with and without IRAC pho¬ 
tometry. The error bars are based on Monte Carlo modeling given 
the photometric errors. There are very few galaxies for which the 
stellar mass is significantly underpredicted when IRAC photome¬ 
try is excluded. It is more common for the no-IRAC fits to biased 
towards higher stellar masses, in at least some cases because of the 
extra contribution to the i^ s -band light by Ho; + [Nil] emission. 
The red points are galaxies whose redshifts place Ha outside of the 
K s band, while the blue points are objects with 2.016 < z < 2.504 
for which significant contamination of the K s measurement by Ha 
emission is possible. 



Without IRAC photometry, we find typical Monte 
Carlo stellar mass uncertainties of 50—60% for CSF mod¬ 
els, which is ~ 1.5 times larger than the typical stellar 
mass uncertainties estimated when including the IRAC 
photometry. When r is allowed to vary, the typical stel¬ 
lar mass uncertainty is 60 — 70%, ~ 2 times larger than 
the typical uncertainty inferred with IRAC photometry. 
Thus, while the lack of IRAC photometry does not ap¬ 
pear to cause a large bias in the stellar mass estimates, 
the inclusion of IRAC photometry significantly reduces 
the formal uncertainties. 

While we have noted the possibility of Ha contamina¬ 
tion of the A' s -band measurements for the ~60% of the 
current sample having redshifts where this is a poten¬ 
tial problem, the total inferred stellar mass for the full 
sample of 72 galaxies differs by only ~ 15% when IRAC 
measurements are included or excluded (in the sense that 
the inferred total mass is higher without the IRAC data 
included). Thus, while the addition of longer-wavelength 
data leads to significantly smaller uncertainties in the es¬ 
timates of the stellar masses of individual galaxies, the 
IRAC data do not appear to be essential for obtaining 
reasonable stellar mass estimates for galaxies at z ~ 2, 
particularly when statistics of large samples are most rel¬ 
evant. 

5.4. Mass-to-Light Ratios versus Wavelength 

It is often asserted that estimates of stellar mass from 
photometric measurements become increasingly reliable 
as one obtains longer wavelength data. We have certainly 
found this to be the case, as discussed above; however, it 


is not true that long-wavelength observations can be used 
to measure stellar masses without understanding the re¬ 
cent star-formation history of the galaxy, even when ob¬ 
servations extend into the rest-frame near-IR as for the 
IRAC observations of z ~ 2 galaxies. This point is illus¬ 
trated in figures El and E3 we have chosen the observed 
I\ s band and the observed IRAC channel 2 (4.5pm) band 
since the former is the longest wavelength measurement 
possible from the ground, and the latter lies close to 
~ 1.5/MR in the rest frame at the mean redshift of the 
current sample (68/72 of the galaxies in the sample have 
been measured in IRAC channel 2), where the peak in the 
SED produced by older stellar populations is expected. 

Figure El shows the absolute magnitude of the galaxies 
based on their observed AT s -band magnitudes, which at 
the median redshift of (z) = 2.29 corresponds to rest- 
frame 0.65pm (i.e., R band). Note that while there is 
clearly a correlation between absolute R-band magni¬ 
tude and stellar mass, there is a large enough range in 
star-formation history among the sample that a measure¬ 
ment of a particular rest-frame luminosity would map 
quite poorly onto stellar mass. The dotted and dashed 
curves show that there is a range in inferred M/L of 
a factor of ~ 70 among the sample. This large fac¬ 
tor applies even for objects that are bright in the I\ s 
band (K s < 20.0), indicated with open circles. The very 
largest stell ar M/L value s approach that of present-day 
galaxies ( Hell et alJl2003lh but the average value is ~ 5 
times smaller, ((M/L)a) = 0.53 (where both M and L 
are in solar units, and L is evaluated in the rest-frame 
R band). The situation is improved at 4.5pm (~ rest- 
frame 1.4pm), as shown in figure ITUl the scatter in the 
stellar M/L is a factor of ~ 15, with indications that the 
scatter decreases for the most luminous galaxies in the 
rest-frame near-IR, which tend to have M/L very close 
to that of present-day galaxies. The scatter in M/L is 
also greatly reduced when the sample is limited to ob¬ 
jects with 7 Z — K s > 3.5 (the reddest quartile), indicated 
with open squares. Still, among the total sample stud¬ 
ied here, the average value is (M/L) i. 4 Mm = 0.38, roughly 
2.5 times smaller than typical in the local universe at the 
same rest wavelengths. The local M/L ratios were de¬ 
rived with an IMF containing fewer low-mass stars than 
are found in the standard 0.1 — 1OOM 0 Salpeter IMF, 
used for the cu rren t work. Adopting the “diet” Salpeter 
IMF used in iBell et ahl (1200311 would lower the average 
2 ~ 2 M/L ratios by a factor of ~ 1.5. 

A simple interpretation of these results is that gaug¬ 
ing the importance of current or recent star formation on 
the galaxy luminosities, even in the rest-frame near-IR, 
depends on the observations of the full UV-to-near-IR 
SED of the galaxy. That is, one cannot go directly from 
near-IR luminosity to stellar mass without knowledge of 
the likely contribution of young stars to that luminosity, 
which can only come from measurements at shorter wave¬ 
lengths. The region of the spectrum that is most sensitive 
to the ratio between the current star formation and the 
integral of past star formation is between the UV and the 
visual, conveniently measured by the observed 1Z — K s 
color in the current sample. This point is illustrated in 
figure M which shows that the inferred stellar mass is 
remarkably well-correlated with 1Z — K S . In fact, the cor¬ 
relation between 1Z — K s color and inferred stellar mass 
is almost as significant as the correlation between stellar 
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Fig. 9. — Inferred stellar mass for the best-fit model, versus 
the absolute magnitude at K s (2.155/mi) Note that the scatter 
in stellar mass/light ratios at a given K s luminosity (rest-frame R 
band at the mean redshift of the sample) can be as large as a factor 
of 70. The 9 objects in the sample with K s < 20 (Vega normalized) 
are indicated with open circles surrounding the points, while the 
reddest quartile of the sample with 1Z — K s > 3.5 are indicated 
with open squares. Dashed and dotted lines bracket the range of 
stell ar M/L in the sample, and are given in solar units, evaluated 
at rest-frame R band. The very largest optical stellar M/L seen 
in z rsj 2 galaxies approaches that of galaxies in the present-day 
universe, while the average value is ~ 5 times smaller. 



mass and using a Spearman test, the signifi¬ 

cance of the correlations are 5.9a and 6.4cr, respectively 
- and more significant than the correlation between stel¬ 
lar mass and rest-frame 0.65/im absolute magnitude (5a). 
One implication of this result is that, in the absence of 
IRAC data, it may be more effective to use 1Z — K s color 
rather then K s flux to select z ~ 2 galaxies with large 
stellar masses. Another implication is that there are few 
galaxies that have large inferred stellar masses at z ~ 2 
without showing evidence fo r long star-formation histo¬ 
ries Ccf.. iShaolev et alll2004 1 that are suggested by red 
rest-frame 2000 6500 A colors. 

Of course, caution is required in assessing the true sig¬ 
nificance of correlations that are based on the model 
results and not on directly measured quantities; such 
is the case for all stellar mass estimates. To some ex¬ 
tent the results in this section simply reflect how the 
observed SED is used to infer stellar mass by the mod¬ 
eling software, as a single smoothly declining or contin¬ 
uous episode of star formation that dominates the emis¬ 
sion at all wavelengths from rest-frame UV through near- 
IR. In such a framework, the ratio of rest-frame UV to 
optical or near-IR light indicates the stellar M/L in a 
fairly robust manner (see section im We can obtain 
a different result for the stellar mass if we instead un¬ 
derstand the galaxy UV-to-near-IR SED as the super¬ 
position of a young, roughly continuous episode of star 
formation and an old burst that peaked sometime in the 
past. In the most extreme case of a maximally-old burst 
and a very young (t s f < 10 Myr) continuous episode, 
the old burst will dominate the rest-frame near-IR light, 
the young episode will dominate the emission at rest- 
frame UV wavelengths, and the maximum possible stel¬ 
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Fig. 10.— Inferred stellar mass for the best-fit model, versus the 
absolute magnitude at 4.5pm, equivalent to the rest-frame 1.4pm 
luminosity at the mean redshift of the sample. As for figured the 
objects with K s < 20.0 are circled and those with 7 Z — K s > 3.5 are 
indicated with open squares. Dashed and dotted lines bracket the 
range of stellar M/L in the sample, and are given in solar units, 
evaluated at rest-frame 1.4pm. While the scatter in stellar M/L 
is very much reduced compared to that measure in the K s band, 
it can still be as large as a factor of 10 at a given 4.5pm lumi¬ 
nosity. The z ~ 2 sample spans a range of a factor of 15, with 
the most massive galaxies having attained a near-IR M/L com¬ 
parable to present-day galaxies, but with typical galaxies having 
(M/L)!, 4 um — 0.4, a factor of ~ 2.5 smaller. There is much less 
scatter in near-IR M/L for objects with R, — K s > 3.5. 



lar mass will be inferred from the rest-frame near-IR lu¬ 
minosities. While this scenario is rather unrealistic, it 
does provide an extreme upper limit on the stellar mass 
implied by the rest-frame near-IR light. Less extreme 
multiple-component models provide a natural represen¬ 
tation of the episodic star-formation histories that arise 
in a hierarchical model of galaxy formation, and will al¬ 
ways result in larger stellar mass estimates than those 
derived from the single-component modeling presented 
thus far in this work. 

In order to set a rough upper limit on the stellar masses 
implied by the broad-band SED s, we use d a technique 
very similar to that of lDaddi et all (' 2004H1 . which is de¬ 
signed to infer the maximum M/L ratio from the rest- 
frame near-IR data points. This method accounts for 
most of the light at those wavelengths with a maximally- 
old burst, on top of which is superimposed a young stellar 
population that accounts for most of the rest-frame UV 
emission. The ages of the maximally old models range 
from 2.2-4.3 Gyr for the galaxies in our sample. We first 
fit a young ( t s f < 10 Myr) CSF model to the observed 
photometry, normalizing this model to the 1Z magnitude. 
We then subtracted this young model from the observed 
photometry, normalized the maximally old r = 100 Myr 
burst to the residual IRAC magnitudes, and summed 
the stellar mass from the two components. This sum is 
dominated by the mass of the old component, and is typ¬ 
ically ~ 3 times larger than the single-component mass. 
The correlation between M, 4 Mm and stellar mass remains 
significant at the same level as for the single-component 
models. While 1Z—K S color is still significantly correlated 
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Fig. 11. — The relationship between inferred stellar mass and 
the observed 1Z — JC color for the sample of 72 galaxies with (z) = 
2.26 ±0.30. There is a very strong correlation between inferred M* 
and 1Z — K s color, particularly for 1Z — K s >3.5. Again, objects 
with K s <20.0 are circled. 
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Fig. 12. — The relationship between the absolute magnitude in 
the observed 7Z band (equivalent to rest-frame 2100 A at the mean 
redshift of the sample) and inferred stellar mass. A Spearman test 
indicates no correlation of these two quantities. There is a similar 
lack of correlation between M* and the absolute magnitude at rest- 
frame ~ 1500 A. Again, objects with K s < 20.0 are circled and 
those with 1Z — K s >3.5 are indicated with open squares. 


with stellar mass (at the ~ 4 a level), the correlation is 
less significant than for the case of the single-component 
models discussed above. Specifically, there is much more 
scatter among objects with blue TZ — K s colors, for which 
the M/L can change the most using this type of model¬ 
ing. For objects with 7 Z — K s > 3.5, the inferred M/L 
and the correlation between stellar mass and 1Z — K s 
color are fairly insensitive to the type of model used. 

It is easy to see that the amount of additional stellar 
mass that could be “hidden” by the current star forma¬ 
tion episode is strongly dependent on the observed ratio 
of UV to optical/near-IR light. For the bluest objects, 
the ratio of true stellar mass to inferred stellar mass could 
be > 5 while redder galaxies (except in extreme cases of 
dust obscuration) could not hide much additional stellar 
mass. For the galaxies in the sample with the largest stel¬ 
lar masses derived from single-component fits, the two 
component modeling provides stellar mass upper limits 
only 10—40% higher than the original estimates. Thus, in 
general, objects which already have large inferred stellar 
masses and old inferred ages cannot be much altered by 
more extreme assumptions, while objects that are com¬ 
pletely dominated by the current episode of star forma¬ 
tion (these would tend to be galaxies with very young 
inferred ages) could be more massive by factors as large 
as ~ 5 — 10. The typical BX/MD galaxy, which is be¬ 
tween these extremes, probably has a true stellar mass 
that is a factor of < 3 larger than the value inferred from 
population synthesis modeling. Finally, this decomposi¬ 
tion of the SED into young and maximally-old compo¬ 
nents provides a fairly robust (but extreme) upper limit 
to the stellar mass consistent with the observed photom¬ 
etry, since more realistic, longer current star-formation 
episodes, would make an increasing contribution to the 
rest-frame near-IR light, and more realistic past bursts 
would occur at z < oo and have smaller M/L ratios. 

5.5. AGN 


As described above, there is one object, MD157, whose 
spectrum shows clear signatures of AGN activity in the 
far-UV (strong C IV and He II emission as well as a broad 
Lyman a line). In addition, MD157 is extremely bright 
in the IRAC-4 band (AB mag of 18.51; see Table 2), such 
that obtaining a reasonable fit to its SED using stellar 
populations required excluding the 8/rm data point. Two 
other objects, MD174 and MD94, have rest-frame opti¬ 
cal spectra indicative of AGN (both have Ha/[Nil] ~ 1, 
and broad Ha with vfwhm > 1000 km s -1 ). All three 
objects exhibit significantly redder than average 7 Z — K s 
colors, ranging from 4.46 — 5.15 and best-fit stellar masses 
ranging from 8 x 10 1O M 0 to 4.5 x 10 11 M S . 

An important question is how much the presence of the 
AGN affects the stellar mass estimates. The large stellar 
masses are inferred because of the bright rest-frame near- 
IR photometry and the very red rest-frame UV/optical 
and UV/near-IR colors, but if a dusty AGN were dom¬ 
inating the rest-IR part of the SED, these stellar mass 
inferences could be erroneous. A large fraction of other 
similarly bright and red z /L 2 galaxies having relevant 
spectroscopic measurements also indicate the presence of 
AGN. For example, 4 of the 7 gal axies w i th spectros copic 
redshifts fr om the FIRES surv ey mx^k SJEHqS) pre¬ 

sented by Ivan Dokkum et all (f2003L l2004fl exhibit UV 
and/or optical signatures of AGN activity (high ioniza¬ 
tion UV lines, broad Ha and Ha/[Nil] line ratios near 
unity), similar to MD157, MD174, and MD94. It is diffi¬ 
cult to know whether AGN activity leads to erroneously 
large stellar mass estimates, or whether it is preferen¬ 
tially objects that have already accumulated large stel¬ 
lar masses that tend to harbor significant AGN activity. 
When the AGN is not obviously energetically dominant 
in the observed-frame mid-IR, we suspect that the stellar 
mass estimates are probably reasonable; for objects like 
MD157, which is obvious as an AGN in both the rest- 
frame far-UV and near-IR, it is probably safest to assume 
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that the stellar mass has been over-estimated using SED 
fitting. 

There are other objects that appear to have fluxes at 
8/xrn that are significantly in excess of that expected from 
the SED modeling (e.g., BX756, BX535, MD69, BX505, 
BX581, BX561, BX523) but where we do not have any 
obvious spectroscopic evidence for the presence of an 
AGN. It is conceivable that these rnid-IR excesses could 
be due to obscured AGN activity; deep Chandra images 
of this field could help resolve this question. 

6. DISCUSSION 

We have seen above that modeling of the broad-band 
SEDs of z ~ 2 star forming galaxies has significant 
limitations, even with spectroscopic redshifts and when 
longer-wavelength data are available for direct measure¬ 
ment of the rest-frame near-IR light. We have shown 
that reasonably consistent stellar mass estimates are ob¬ 
tained independent of the exact form of the star forma¬ 
tion history, but that this consistency should be viewed 
with caution because of the fact that a robust episode of 
star formation can change the inferred M/L by a factor 
of several, even at rest-frame ~ 1.5 /im. One obtains in¬ 
creasingly reliable estimates of the total stellar mass of a 
galaxy as the integral of star formation over the lifetime 
of the galaxy exceeds that in the past several hundred 
Myr by a factor of a few. In the most reliable cases, 
where the UV light is far surpassed by the emission at 
rest-frame optical and near-IR wavelengths, the inferred 
stellar mass does not depend very strongly on whether 
a smooth or episodic prescription is used to model the 
way in which the integrated stellar mass accumulated as 
a function of time. When the rest-frame UV luminos¬ 
ity is more comparable to rest-frame optical and near- 
IR luminosities, however, it is inevitable that the stellar 
masses inferred from photometric techniques will reflect 
the typical stellar mass produced in an episode of star for¬ 
mation, perhaps caused by a merger or accretion event 
over the course of the last few hundred Myr, and not the 
total stellar mass. A corollary to the above arguments 
is that it is most likely one would measure large stel¬ 
lar masses (using the type of modeling described in this 
paper) when the current rate of star formation is signifi¬ 
cantly smaller than the past average. Otherwise, most of 
the rest-frame near-IR and optical light will be ascribed 
to the current star-formation episode, which has much 
smaller M/L even in the IRAC bands. These are not new 
argum ents (see, e.g.. lPanovich et aTbOOltfShaplev et all 
Eiomlh though the inclusion of IRAC wavelengths to con¬ 
strain the rest-frame near-IR emission places them on 
a more robust footing, and the larger range of galaxy 
properties observed in the current sample highlights the 
important trends more clearly. 

The concept of episodic star formation informs the 
above discussion. While such variable star-formation his¬ 
tories appear a reasonable representation of the build¬ 
up of stellar mass in galaxies at high redshift, sim¬ 
ply based on general consideration of the hierarchical 
model of galaxy formation in which merging (and merger- 
induced star formation) plays a large role, we also find 
direct evidence from the galaxies in our sample. Specif¬ 
ically, we can obtain interesting constraints on the red- 
shift at which the onset of star formation occurred for 
the galaxies in our sample, despite large systematic un¬ 


certainties in the stellar population ages (section 14.311 . 
The maximum stellar population age is obtained for a 
given set of rest-frame UV to near-IR colors when a con¬ 
tinuous star-formation history is used, so we can use 
the CSF results to infer the earliest redshifts at which 
the c ur rent episodes of star formation could begin (cf., 
iFerguson et al.ll200ilh The median CSF age for the sam¬ 
ple is ~ 600 Myr. Using the redshift and CSF age distri¬ 
butions we find that the current episode of star formation 
had not begun yet at z = 3 for at least ~ 65% (46/72) of 
the galaxies in the sample. The true fraction is probably 
higher, since CSF models provide only an upper limit 
on the age of the stellar population dominating the cur¬ 
rent star-formation episode. Therefore, we find that the 
majority of galaxies in our sample have young ages for 
the current episode of star formation compared to the 
age of the Universe at z ~ 2, allowing the possibility of 
much older previous bursts of star formation. Relatively 
young inferred ages would naturally arise in a picture of 
episodic star-formation. 

A related point is that selecting on the basis of far- 
UV properties, as we have done for the BX/MD sample, 
connects only obliquely to stellar mass. One might ex¬ 
pect that galaxies would have to attain a particular total 
mass threshold before they could support star formation 
detectable at the apparent magnitude limit of the sur¬ 
vey, but once that threshold is exceeded, a combination 
of extinction (whic h is pro bably sta r -form ation rate de¬ 
pendent; e.g.. lAdelberger fc Steidell ( 2000 )1 and scatter 
in SFR per unit total (or stellar) mass might “random¬ 
ize” the mix of galaxy properties included in the sample. 
Indeed, as shown in figure El we find that there is no 
hint of a correlation between Mjjy and inferred stellar 
m ass, as was fou nd pre viously in the z ~ 3 LBG sample 
of lShanlev et alJ ll 20011b 

These issues make it particularly difficult to measure 
the total stellar mass that has formed by a particular 
redshift llDickins on et al.ll2003t iGlazebrook et' all 12004 
iFontana et al J12Q03|), or to make a detailed comparison 
of the stellar mass of even similarly-selected galaxies at 
different redshifts. For example, analysis of the spa¬ 
tial clustering of z ~ 3 LBGs and z ~ 2 BX galax¬ 
ies l|Adelberger et, alJ 1 2004 all suggests that correlation 
length of LBG descendants at z ~ 2 would be similar 
to the observed BX correlation length, with the impli¬ 
cation that z ~ 2 LBG descendants reside in dark mat¬ 
ter halos in the same mass range as those hosting BX 
galaxies. It is therefore of interest also to compare the 
stellar mass distribution of spectroscopically-confirmed 
BX/ MD o b jects w ith that of z ~ 3 LBGs from the sam¬ 
ple of lShaolev et all (1 2001 1. A direct evolutionary com¬ 
parison between the stellar mass distributions for z ~ 2 
and z ~ 3 UV-selected samples is not straightforward 
for many reasons. Most importantly, while there is little 
doubt that the objects at z ~ 2 and z ~ 3 with large 
inferre d ma sses and ages really are massive and old (see 
section IQ) . the modeling results likely represent lower 
limits on the stellar masses for more typical objects in 
both samples, where the current star formation is pro¬ 
ducing enough optical and near-IR light to potentially 
mask previous generations of stars. If star formation is 
episodic, as seems quite likely on many grounds, then it 
seems possible in principle that the same galaxy could be 
measured at different times, but the inferred stellar mass 

























20 


might not even be a monotonically increasing function of 
time. Similarly, it might be difficult to trace the evolu¬ 
tion of stellar mass during the epoch when most galaxies 
are still very actively star-forming. 

Bearing all of these concerns in mind, we now summa¬ 
rize the inferred properties of typical UV-selected galax¬ 
ies at z ~ 2. We then turn in some detail to the most 
massive and oldest galaxies in the sample, for which we 
obtain robust total stellar mass measurements, and show 
how these massive galaxies provide plausible progenitors 
for passive galaxies identified at z ~ 1.5 —2.0. Finally, we 
make a general comment about the relationships among 
galaxy populations selected with different photometric 
criteria. 


6.1. The Properties of the “Typical” UV-Selected 
Galaxy at z ~ 2 


The results of this paper suggest that the typical 
0 ~ 2.2 “BX” galaxy has a stellar mass of ~ 3 x 10 10 
M 0 , a number that could possibly be higher by a factor 
of ~ 2 — 3 without affecting the observed SEDs. Al¬ 
though the length of time stars have been forming in in¬ 
dividual ^alaxiesjs not well-determined by the modeling, 
iSteidel et al.l 12005 1 show that galaxies within the signif¬ 
icant over-density at z = 2.30 are older and contain more 
stellar mass on average than those in lower-density en¬ 
vironments. It is very clear from these results that even 
at z ~ 2 the star-formation history is strongly depen¬ 
dent on large scale environment- more so than on the 
particular redshift at which the galaxies are observed. 
The typical current sta r-formation rate is ~ 50 M 0 yr _1 
llReddv & SteideH 12004 *1. although unobscured star for¬ 
mation down to ~ 3 M 0 yr _1 could make it into the 
sample, and we know that among the BX sample lurk 
sub-mm s ources that have impli ed SFR of up to ^ 1000 
M 0 yr _1 llChaoman et all 120041 . UV-selected galaxies 
that are also optically bright ( K s ;$ 20) h ave typical 
inferre d SFR higher by a factor of ~ 3 llReddv et al.l 
l2005al) . 


The space density of galaxies similar to those dis- 
cussed i n this paper is ~ 3 x 10~ 3 Mpc -3 at z zz 2.2 
llReddv et al.ll2005bl) . The co-moving correlation length 
of the U V-selected galaxies at z ~ 2 is r Q = 4.2 ± 0.4 
h _1 Mpc (lAdelberger et al.ll2005l) . implying that the typ¬ 
ical galaxy fives in a dark matter halo of mass ~ 2 x 10 12 
M 0 . Assuming that the galaxies have a baryon frac¬ 
tion that is the same as the universal value, the es¬ 
timated mass in baryons associated with the halo is 
~3x 10 11 M 0 , suggesting that on average roughly 10% 
of the baryons have been converted into stars. Con¬ 
sideration of the expected evolution of the clustering 
properties of the UV-selected galaxies suggests that by 
z ~ 1 they will be found mainly among early-type galax¬ 
ies, and will have completed most of their star forma¬ 
tion. The most massive examples will have completed 
their star formation even earlier, and exhibit significantly 
stronger clusterin g than the average UV-selected galax¬ 
ies at z ~ 2 ijAdelberger et al.l 2 01)51 :. As discussed fur¬ 
ther in section IQ and other recent work, the ~ 2.5 
Gyr between z ~ 2 and z ~ 1 marks the transition 
of a large fraction of massive gal axies f rom strongly 


Shanlev et al.12004 

Cimatti et al. 2004 McCarthy et al. 

20041 iFontana et al. 
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6.2. Massive Galaxies in Rest-UV-Selected, Samples 

The most massive galaxies in the current sample, 
with stellar masses > 10 11 -M 0 , cannot contain signifi¬ 
cant additional stellar mass that is masked by the cur¬ 
rent episode of star formation. We obtain roughly the 
same stellar mass estimate regardless of whether a single¬ 
component model is used to fit the colors, or a young 
current star-formation episode super-imposed on an old 
past burst. In the context of simple, single-component 
models, these galaxies also have SEDs that provide sig¬ 
nificant constraints on the star-formation history. Unlike 
the majority of galaxies in the HS1700 sample, for which 
we find that all or most r values are consistent at the 
95.4% level with the data, we can rule out models with 
t < 200 Myr and r = oo (i.e., CSF) for all but one of 
the galaxies with M* > 10 11 M 0 . More discrimination 
among r values leads to interesting constraints on the 
duration of the star-formation episodes in these galaxies. 
The minimum r = 200 Myr for these five galaxies im¬ 
plies minimum ages of t s f > 600 — 900 Myr. The most 
likely value of r however is 1 — 2 Gyr, based on Monte 
Carlo simulations in which r was allowed to vary. For 
r = 1 — 2 Gyr, the ages inferred for these galaxies range 
from 2 — 3 Gyr, comparable to the age of the Universe at 
0 ~ 2.3. The stellar masses for this subsample of galaxies 
are all well-constrained (< 20% uncertainties), including 
the systematic uncertainty in star-formation history 8 . 

The superior constraint on star-formation history for 
these objects results from two factors: first, the most 
massive galaxies are among the brightest in the K s and 
IRAC bands 9 , resulting in smaller than average TZ — K s 
and IRAC photometric uncertainties. Four of the five 
massive galaxies with significant constraints on r are also 
detected in all four IRAC bands, providing better con¬ 
straints on the shape of the SED in the region dominated 
by longer-lived stars. Secondly, all but one of the objects 
with M* > 10 n M 0 have significantly redder than aver¬ 
age TZ — K s (and 1Z — IRAC) colors 10 , but average or 
bluer than average G — TZ colors. The combination of 
blue to average G — TZ and red 7 Z — K s colors indicates 
that the rest-frame optical and near-IR light cannot be 
ascribed to the massive stars tracing the instantaneous 
star-formation rate and dominating the rest-frame UV 
emission, and therefore must be produced by less mas¬ 
sive stars formed primarily in the past. Since the inferred 
extinction correction will be small when the UV color is 
blue, either the star formation was much higher in the 
past, or the current SFR has been sustained for a very 
long time. If the implied CSF time is older than the age 
of the universe, then r models will be required to fit the 
data well. 

There are other objects in the sample with comparably 
small TZ — K s and IRAC uncertainties, but with weaker 

8 We note that the massive galax ies i n the current sam ple are 
quite similar to those presented by ISliaolev el al. ( 2004 *1, which 
were all drawn from another of our survey fields, Q1623+26. 

9 We note, for comparison with the K20 survey (2imatti_e^alJ 
l200l that only half of the galaxies in our sample with M* > 
10 11 Mq also have K s < 20.0. A deeper apparent magnitude limit 
than K s = 20.0 is evidently necessa ry to g ain a comp l ete ce nsus of 
the most massive galaxies at z ~ 2 IGlazebrook et aT]l2004^ . 

10 Most of the galaxies with M* > 10 11 Mq have best-fit models 
with J — K s > 2.3, implying that they would be selected by the 
FIRES survey. 
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constraints on r. An example is MD94, which has very 
red 1Z — K s and G — TZ colors and is inferred to have a 
stellar mass M > lO n M 0 . In this case, the colors can 
be adequately fit by models spanning the whole range 
of t, including CSF, mainly because the inferred current 
SFR (derived from the extremely red G—7Z color) is high 
enough to have produced the near-IR light in much less 
than the age of the universe, even for CSF. The best- 
fit r model, though, has an age that is a factor of ~ 3 
shorter than the CSF model, and a current SFR 3 times 
smaller (however, the implied initial SFR given the age 
and r values would have been 25 times larger, or 2500 
M 0 yr -1 !). 

We return now to the massive galaxies with superior 
constraints on the single-component star-formation his¬ 
tories. Of course, a smoothly-declining star-formation 
history is not a very likely scenario for the build-up of 
stellar mass, since the star-formation rate is expected to 
be modulated by episodic merger and accretion events. It 
is not really possible to establish whether or not star for¬ 
mation has been sustained continuously at a detectable 
level over a long period of time, even for galaxies whose 
star formation must have begun at very high redshift. 
For example, in general, galaxies for which the r models 
are required for a good fit to single-component models 
can also be adequately fit by assuming a single, early 
burst of star formation to produce the optical/near-IR 
light and a small amount of recent star formation (with 
an associated stellar mass of only ~ 1 — 2% of the early 
burst) that accounts for the rest-frame far-UV light. For 
such two-component models, there is again no unique so¬ 
lution, but reasonable fits to the data can be found for 
models where a large burst with r = lOOMyr occurred at 
z > 4 and the recent (smaller) burst is less than 100 Myr 
old 11 . While we cannot distinguish between smooth or 
episodic star-formation histories, the estimates of total 
stellar mass derived from two-component models agree 
very well with those obtained from the simple single¬ 
component fits. For these massive and old galaxies, the 
inferred total stellar mass is therefore a robust quantity. 

In an analysis of the stellar populations of UV-selected 
2.0 < z < 3.5 galaxies in the Hubble Deep Field (HDF- 
N) iPanovich et all : 200 11 found no comparably massive 
and old objects. The small ratio between the typical 
inferred star-formation timescale for these HDF LBGs 
and the time interval between 2.0 < z < 3.5 was used 
to motivate a scenario of multiple, intermittent star- 
formation episodes with a duty cycle of < 1 Gyr. As 
described in the beginning of section M most of the cur¬ 
rent sample can in fact be described in this manner. 
In retrospect, i t is n ot surprising that the sample of 
IPanovich et, all (I2001H did not contain any M* > 10 11 
M 0 galaxies in the 5.0 arcmin 2 HDF-N field, despite 
a similar optical magnitude limit in the HDF spectro¬ 
scopic sub-sample. We found only 6 with spectroscopic 
redshifts in the entire 65 arcmin 2 WIRC field. Since the 
spectroscopic completeness of our sample is ~ 22%, the 
approximate surface density of objects in the redshift 
interval z ~ 2.3 ± 0.4 is ~ 0.5 arcmin -2 , so that one 
would expect < 2 such massive and old galaxies in an 

11 We note that such models would imply SFRs that averaged 
~ 1500 Mq yr -1 over the first 200 Myr of the galaxy’s life, e.g. 
between 2 = 5 and z = 4.4. 


area with the size of the H DF-N giv en spectroscopic in¬ 
completeness. In fact, lLabbe et all ( 20031 ) have pointed 
out the discrepancy in the number of galaxies with red 
near-IR colors in the HDF-N versus their FIRES field in 
the HDF-S, with the HDF-N having only 10-20% of the 
numbers in the HDF-S. Given that the galaxies with la rge 
stellar mass are s t rongly clustered llDaddi et alJ l2f)04bl) . 
lAdelbererer et aD 1201)5 11 the variance between small sur¬ 
vey areas is expected to be large for such rare objects. 

6.3. UV-Selected Progenitors of Massive and Passive 
Galaxies? 

While massive and mature galaxies represent a small 
fraction by number of the BX/MD population (6/72 with 
M > lO n M 0 , 5 of which also have tsF ~ the age of the 
Universe), their estimated space density, ~ 10 -4 Mpc -3 , 
is similar to the space density of the massive star-f orming 
galaxie s from the K20 survey presented bv lDaddi et all 
( 2004b ) at 1.7 < z < 2.3, and also the massive galax¬ 
ies with old and passiv e ste llar populations presen ted 
by IGimatti et alJ |2004) and IMcGarthv et alJ 120041) at 
1.3 < z < 2.0. In both of the latter papers, the authors 
suggest that the progenitors of the old and passive galax¬ 
ies are to be found among submillimeter- and near-IR- 
selected samples at z > 2. Here we show that the most 
massive and mature objects in the 2.0 < z < 2.6 UV- 
selected sample could quite plausibly be the progenitors 
of massive and passive systems at 1.3 < z <2.0. To make 
this comparison, we adopt star formation histories that 
are consistent both with those used to fit the near-I R 
selecte d m assive and passiv e galax ies by IGimatti et alJ 
(2004) and iMcCarthv et~ahl 120041) and with the SEDs 
of all 5 of the BX/MD massive and old galaxies at 
0 ~ 2.3. We then evolve the se models forw ard in time 
to the mean redshifts of the IGimatti et all 120041) and 
IMcGarthv et all 120041) samples, compute the average K s 
magnitudes, R — K s and I — K s colors, stellar masses, 
and star-formation rates, and compare them with the 
prop erties o f the passi ve galaxies. 

IGimatti et al . ( 120041 ) present four galaxies at (z) = 1.7, 
with (K s ) = 18.7, (R — K s ) = 6.3, stellar masses rang¬ 
ing from 1 — 3 x 1O 11 M 0 , and a space density of ~ 
10 -4 Mpc -3 . When a declining model with r = 0.3 Gyr 
is used to fit the colors, the associated ages and forma¬ 
tion redshifts are ~ 2 Gyr, and z ~ 4, respectively. If 
we apply a r = 0.3Gyr model to the massive BX/MD 
SEDs at (z) = 2.3, and then evolve them forward to 
z = 1.70, they would have (K s ) = 19.4, (R — K s ) = 6.2, 
(M*) = 2.2 x lO n M 0 , ages ranging from 1.9-2.3 Gyr, 
and SFRs of < 1 M 0 yr , making them fainter than 
TZ = 25.0. We find a similar level of agreement between 
the BX/MD models e volved to z = 1.5 and th e parame¬ 
ters of the 5 objects in IMcGarthv et ah] 1 2004 ) (their ta¬ 
ble 1) for which a r = 0.5 Gyr model was used to fit the 
spectra and colors. Thus, extrapolating star formation 
histories for the massive BX/MDs forward in time would 
make them plausibly consistent with the passive galax¬ 
ies in the near-IR selected samples observed at somewhat 
lower redshift. Equivalently, if the BX/MD examples are 
interpreted as massive old bursts with a small amount of 
recent star formation superposed, by z ~ 1.7 they would 
have essentially the same characteristics as the passive 
galaxies provided that the current episode of star forma- 
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tion (i.e., as observed at z ~ 2.3) lasts for less than the 
~ 900 Myr between 2 ~ 2.3 and z ~ 1.7. 

The specific values of the timescale over which star for¬ 
mation occurs, and the redshift at which star formation 
begins, are not well-determined because of the system¬ 
atic uncertainties inherent in stellar population synthesis 
modeling. However, we have demonstrated that, for rea¬ 
sonable assumptions (or at least consistent with those 
used by other workers), the most massive and evolved 
galaxies in the UV-selected sample are likely to evolve 
into largely passive systems by z ~ 1.7, with masses, 
ages, inferred formation redshifts, colors, and space den¬ 
sities similar to those of the massive and old stellar pop¬ 
ulations presented in the K20 and GDDS surveys. While 
additional progenitors of the 1.3 < z < 2.0 massive and 
passive systems may be found among submillimeter- and 
J — K color-selected samples, the UV-selected sample at 
z ~ 2.3 must contain a significant fraction of the objects 
t hat will look like e l liptical galaxies at 1.3 < z < 2.0. 

lAdelberger et al. t 2004a ) have reached similar conclu¬ 
sions about the fate of rest-UV-selected galaxies, using 
independent lines of argument. On the basis of a com¬ 
parison of the observed clustering properties and galaxy 
number density at z ~ 2 with dark matter halos in nu¬ 
merical simulations, Adelberger et al. argue that the BX 
galaxies are consistent with typical total halo masses of 
^ 10 12 Mq and that by z ~ 1 they will have a clustering 
strength similar to that of galaxies with early type spec¬ 
tra in the DEEP2 redshift survey results of iCoil et al.l 
(20041. In this context, the most massive BX galaxies 
at z ~ 2 would be among the more extreme (in terms 
of cl ustering prop erties ) early type galaxies by z ~ 1. 
lAdelberger et al.1 (1200 57) has also shown that optically 
bright ( K s < 20.5) rest-UV-selected galaxies at z ~ 2, of 
which the most massive galaxies in the present sample are 
a subset, are as strongly c lustered as their rest-opti cally 
selected counterparts Cch. lDaddi et alJl2003 . 2004b ). 

6.4. Selection Biases and “Populations” 

We consider it likely that real galaxies actually evolve 
in and out of high redshift galaxy samples selected in 
various ways. For example, the massive galaxies in the 
BX/MD sample could easily have been sub-mm galaxies 
at z ~ 4 (when some versions of their star-formation his¬ 
tories predict high enough SFRs for them to have been 
easily detectable at 850 /im, although most would prob¬ 
ably not have been /iJy radio sources because their red¬ 
shift is too high), and may have oscillated between being 
quiescent and being UV-selectable depending on whether 
or not there was a trace of relatively unobscured star for¬ 
mation superposed on the more mature stellar popula¬ 
tions. 

As discussed by Sha plev et alJ 120041) and 
lAdelberger et alJ <l2004a . 120051) . the most massive 
and optically brightest galaxies in the UV-selected 
samples appear to have similar masses, chemical 
abundances, star-formation rates, inferred ages, and 
clustering properties as galaxies selected in the near-IR 
by either apparent brightness or optical/IR color. This 
comparison shows that a UV-selected survey is capable 
of finding examples of galaxies that are very similar to 
those selected in other ways; it is not a demonstration 
that a UV-selected survey “finds” the vast majority of 
galaxies at z ~ 2. That rest-UV selection is incomplete 


with respect to all galaxies present at similar redshifts 
has been shown clearly by the near-IR selected GDDS, 
K20, and FIRES surveys as discussed above, as well 
as by sub-mm selected surveys (e.g., iGhanman et alJ 
(20031). The point is that it is relatively easy for 
objects to come and go from a UV-selected surveys 
on timescales much shorter than a Hubble time at the 
relevant redshifts, independent of their past (or future) 
star-formation histories. Most of the e xisting evidence, 
and s ome theoretical studies (e.g., iNagamine et alJ 
lt200l ). suggests that at z ~ 2, relatively massive 
galaxies (M* x) 10 10 M 0 ) are “UV-selectable” using 
criteria similar to those employed in the current sample 
(7, 50% of the time. 

While such samples would be significantly incomplete 
with respect to all galaxies at a given stellar mass, they 
may not be entirely missing any type of galaxy. The 
fraction that happens to have enough UV light to be se¬ 
lected has significant advantages: spectroscopy is easy 
in both the far-UV and optical (rest frame). These ob¬ 
servations allow for more detailed astrophysical studies, 
and an evaluation of the dependen ce of galaxy proper- 
ties on large-scale environment as in lSteidel et al.l : 2005 '. 
The recent success of t he spectroscopic follow -up of sub- 
mm-selected galaxies IGhanman et, alJ 120041) has bene¬ 
fited from exactly the same advantage: even UV light 
that is inconsequential compared to the bolometric lu¬ 
minosity allows spectroscopic measurements, leading to 
more detailed astrophysical insight. 

7. SUMMARY 

We have presented a detailed analysis of the broad¬ 
band SEDs of a sample of 72 rest-UV selected, spec¬ 
troscopically confirmed, star forming galaxies at z = 
2.3 ± 0.3 in a ~ 8.5 by 8.5 arcmin field centered on 
HS1700+643. This sample of galaxies is broadly rep¬ 
resentative of galaxies selected with well-calibrated pho¬ 
tometric pre-selection using rest-UV colors: all of the 
spectroscopic data were obtained without consideration 
of the longer-wavelength properties. We used extensive 
ground-based photometry from 0.35-2.15 fin r and deep 
Spitzer/IRAC data probing the rest-frame near-IR, to¬ 
gether with information from optical and near-IR spec¬ 
troscopy, to produce an unprecedented data set that al¬ 
lows the most robust statements to date concerning the 
physical state of galaxies at z ~ 2, in particular their 
stellar masses. 

Our main conclusions include the following: 

1. Using population synthesis to model the stellar pop¬ 
ulations, we find that the inferred stellar mass distribu¬ 
tion of the sample is reasonably well-described by a log¬ 
normal distribution: (logjj^) = 10.32 ± 0.51 and that 
the result is essentially unchanged if a constant star for¬ 
mation history is assumed instead of the best-fitting ex¬ 
ponentially declining star-formation history. The typical 
uncertainty in stellar mass, including variations in the as¬ 
sumed form of the star-formation history, is ~ 30 — 35%. 
Sub-samples of the galaxies with brighter observed K s 
(2.15/im) flux have larger inferred stellar mass; 8% of the 
sample has inferred stellar mass > 10 11 Mq. We found 
that fitting SEDs to simultaneously determine photomet¬ 
ric redshifts and stellar population parameters leads to 
large uncertainties in both, and tends to over-estimate 
the inferred stellar mass; thus we use only the sample 
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with spectroscopic identifications. 

2. The use of the Spitzer /IRAC fluxes, which mea¬ 
sure the rest-frame near-IR light from the galaxies and 
thus should yield the strongest available constraints on 
stellar mass, provides a marked improvement in the esti¬ 
mates of the stellar mass of individual galaxies, reducing 
the formal uncertainties by a factor of 1.5 — 2 relative to 
estimates based on photometry that extends only to the 
ground-based K s band. However, we find that the longer 
wavelength data do not change the results significantly: 
the total stellar mass inferred for the entire sample of 72 
galaxies is changed by only 15% depending on whether 
the Spitzer data are used for the modeling or not. At the 
redshifts spanned by our sample, there is a tendency to 
over-estimate slightly the stellar mass when the K s band 
is the longest wavelength observed, because of contami¬ 
nation in the K s band from Ha emission. 

3. The lack of strong constraints on the form of the 
star-formation history remains a problem for most of the 
galaxies in the sample, even when using simple, single¬ 
component models. Ages, exponential time constants, 
and extinction/SFR are largely degenerate with respect 
to each other; however, as found in previous work, the in¬ 
ferred stellar masses remain largely unaffected. External 
constraints can be used to rule out parts of parameter 
space that are otherwise allowed by the SED fits. 

4. While the inferred stellar mass is correlated with 
the flux at 2 /on, the inferred stellar mass-to-light ratios 
can vary by a factor of ~ 70 at a given AT s -band lumi¬ 
nosity (rest-frame R at the mean redshift of the sample). 
Even at an observed wavelength of 4.5/rm (rest-frame 1.4 
/rm) the inferred M/L varies by a factor of 15. At both 
wavelengths (rest-frame 0.65 and 1.4 /mr), the highest 
inferred M/L are similar to those of local galaxies, but 
the average M/L is a factor of at least 5 and 2.5 smaller, 
respectively. Thus, at these redshifts, even rest near-IR 
luminosities are by themselves a relatively poor proxy for 
stellar mass; the entire SED, from the UV to the near-IR, 
is essential. In particular, the conversion from rest-frame 
near-IR luminosity to stellar mass is strongly correlated 
with the observed 1Z—K S color (rest-frame 2000—6500 A) 
which in the context of the models considered is a crude 
measure of the ratio of past to current star formation. 

5. Even with IRAC rest-frame near-IR information, 
there remains the possibility of significant systematic un¬ 
derestimates of the stellar mass using single-component 
SED-fitting techniques; the magnitude of this effect de¬ 
pends sensitively on the shape of the SED. For objects 
that are red in their UV/optical and UV/near-IR col¬ 
ors, the mass inferred from SED fitting is unlikely to 
underestimate the stellar mass by more than a few tens 
of percent, while the bluest objects in the sample could 
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have total stellar masses that are as much as a factor 
of ~ 5 — 10 larger than inferred (in the case of the most 
extreme assumed star-formation histories). A typical ob¬ 
ject in our sample could be “hiding” up to a factor of ~ 3 
more stellar mass than inferred by the modeling. Thus, 
for star-forming galaxies at high redshift, reasonably pre¬ 
cise estimates of stellar mass are possible only for objects 
that are nearly finished forming stars, or caught in a rel¬ 
atively quiescent state between episodes of star forma¬ 
tion. All other estimates should probably be considered 
as lower limits, possibly reflecting the typical stellar mass 
formed during a star-formation episode rather than the 
integrated stellar mass. 

6. We have shown that the most massive galaxies in 
the UV-selected sample, with M* > 10 11 M 0 , could plau¬ 
sibly evolve by z ~ 1.7 into the high redshift “passive” 
galaxies recently identified from near-IR selected spec¬ 
troscopic surveys. The most likely star-formation histo¬ 
ries for these UV-selected examples indicate that they 
began forming stars at very early times and that their 
star-formation rates were considerably higher in the past, 
quite similar to what is found for the near-IR selected 
galaxies. Taken together with the similarities in space 
density, stellar mass, chemical abundances, and cluster¬ 
ing properties, galaxies with similarly large stellar masses 
should probably be treated as a single “population” re¬ 
gardless of the method used to discover them. In most 
cases, the UV light which brings them into the sample 
represents a small fraction of the stellar mass assembly. 
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TABLE 1 

HS 1700+64 Imaging Observations 


Band 

Aeff(Mm) 

Telescope a 

Field Size 

Integration 

Seeing/Image Size 

Photometric Depth b 

U n 

0.36 

WHT/Keck I 2001 May 

15'. 3 x 15'. 3 

23400/3600 

0'.' 89 

26.4 

G 

0.47 

WHT 2001 May 

15'. 3 x 15'. 3 

9900 

0'.' 83 

26.5 

TZ 

0.68 

WHT 2001 May 

15'. 3 x 15'. 3 

7800 

0778 

25.7 

K s 

2.15 

P200 2003 June/Oct 

8'. 5 x 8'. 5 

39720 

0784 

24.0 

IRAC-1 

3.6 

Spitzer 2003 Oct 

5’ xlO’ 

39600 c 

2702 

24.6 

IRAC-2 

4.5 

Spitzer 2003 Oct 

5’ xlO’ 

39600 c 

2738 

24.3 

IRAC-3 

5.8 

Spitzer 2003 Oct 

5’ xlO’ 

39600 c 

2736 

23.0 

IRAC-4 

8.0 

Spitzer 2003 Oct 

5’ xlO’ 

39600 c 

2756 

22.8 


a WHT: William Herschel 4.2m telescope with Prime Focus Imager; Keck I: Low Resolution Imaging Spectrometer, blue 
side. P200: Palomar 5.1m telescope with Wide Field Infrared Camera; Spitzer: Spitzer Space Telescope/Infrared Array 
Camera 

b AB magnitude for a 5cr detection in a 2" diameter aperture. For IRAC bands, the magnitude for a 5cr detection using 
PSF-fitting and local background estimation as described in the text 
c Total exposure time in the deepest part of the mosaic. In general, the deepest regions of the IRAC-1 and IRAC-3 
bands are in different parts of the field from those deepest in IRAC-2 and IRAC-4. 



TABLE 2 

HS 1700+64 Galaxies with Spectroscopic Redshifts 1 


to 

o* 


Name 

o( J2000) 

<5( J2000) 

z 

77 ab 

(G - 77.) ab 

(f/„ - G)ab 

77 - K s b 

(IRAC-I)ab 

(IRAC-2)ab 

(IRAC-3)ab 

(IRAC-4) A b 

BX627 

17:00:34.435 

64:11:29.699 

1.478 

24.44±0.14 

-0.104=0.05 

0.124=0.05 

2.364=0.51 


23.634=0.18 



BX738 

17:01:16.150 

64:12:36.684 

1.710 

23.934=0.12 

0.014=0.04 

0.494=0.09 

3.294=0.15 

22.114=0.10 

22.154=0.10 

22.564=0.29 

23.074=0.49 

BX756 

17:00:55.146 

64:12:54.403 

1.738 

23.21±0.11 

0.33±0.05 

0.58±0.09 

2.414=0.23 


21.824=0.10 


21.554=0.15 

BX681 

17:01:33.762 

64:12:04.278 

1.740 

22.04±0.05 

0.194=0.01 

0.40±0.02 

2.994=0.08 

21.024=0.10 

21.084=0.10 

21.164=0.18 

21.354=0.14 

BX843 

17:01:15.627 

64:13:52.891 

1.771 

24.614=0.16 

0.184=0.08 

0.514=0.11 

3.264=0.31 

22.954=0.34 

22.804=0.10 



BX470 

17:01:26.563 

64:09:31.993 

1.841 

24.26±0.14 

0.38±0.07 

0.75±0.13 

4.204=0.11 

21.444=0.10 

21.454=0.10 

21.304=0.10 

21.694=0.13 

BX846 

17:01:02.925 

64:13:54.470 

1.843 

23.914=0.12 

0.094=0.04 

0.484=0.09 

3.184=0.18 


22.134=0.10 


22.154=0.33 

BX574 

17:00:46.904 

64:10:49.536 

1.846 

24.10±0.14 

0.244=0.07 

0.734=0.10 

3.434=0.19 

22.314=0.10 

22.274=0.10 

22.064=0.22 

22.994=0.36 

BX1043 

17:01:17.328 

64:15:29.170 

1.847 

23.54±0.12 

-0.014:0.05 

0.53±0.09 

2.344=0.42 


22.784=0.24 



BX1087 

17:00:59.047 

64:15:55.195 

1.871 

23.10±0.11 

0.124=0.04 

0.594=0.08 

2.854=0.15 


21.934=0.10 



BX1032 

17:00:52.827 

64:15:23.792 

1.883 

23.764=0.12 

0.184=0.04 

0.624=0.06 

3.044=0.18 


22.374=0.10 


22.594=0.30 

BX782 

17:01:11.146 

64:13:14.360 

1.930 

23.94±0.12 

0.05±0.04 

0.59±0.09 

2.974=0.22 

22.184=0.15 

22.354=0.11 



BX530 

17:00:36.859 

64:10:17.376 

1.942 

23.05±0.11 

0.214=0.05 

0.694=0.05 

2.804=0.15 

21.734=0.10 

21.884=0.10 


22.194=0.33 

BX807 

17:00:54.909 

64:13:28.593 

1.971 

24.834=0.16 

0.08±0.08 

0.454=0.11 

3.044=0.36 


23.494=0.10 



BX903 

17:00:54.431 

64:14:14.761 

1.976 

24.394=0.14 

0.184=0.06 

0.504=0.11 

2.994=0.28 


22.804=0.10 



BX536 

17:01:08.939 

64:10:24.949 

1.977 

23.00±0.11 

0.21±0.05 

0.79±0.05 

3.274=0.10 

21.234=0.10 

21.204=0.10 

21.234=0.10 

21.834=0.12 

BX412 

17:01:20.109 

64:08:37.363 

1.981 

23.584=0.12 

0.314=0.06 

0.634=0.06 

3.314=0.15 

22.134=0.10 

22.234=0.10 

22.324=0.14 


BX526 

17:00:58.040 

64:10:15.947 

2.018 

23.99±0.12 

0.00±0.04 

0.45±0.09 

2.254=0.42 

22.90±0.24 

22.944=0.23 

23.054=0.34 


BX1075 

17:00:51.215 

64:15:47.515 

2.061 

24.83±0.16 

0.30±0.09 

0.89±0.18 

2.924=0.54 


24.464=0.48 



BX625 

17:01:11.493 

64:11:29.242 

2.078 

24.524=0.16 

-0.034=0.07 

0.49±0.06 

2.964=0.54 

23.284=0.10 

23.284=0.10 



BX592 

17:01:02.703 

64:11:03.422 

2.099 

24.874=0.16 

0.384=0.09 

1.114=0.19 

3.014=0.36 

23.354=0.11 

23.434=0.12 



BX535 

17:01:07.597 

64:10:26.403 

2.113 

25.164=0.16 

0.574=0.12 

1.554=0.37 

3.454=0.37 

23.04±0.14 

22.894=0.17 

22.824=0.22 

22.784=0.33 

BX1012 

17:00:47.578 

64:15:11.701 

2.115 

23.72±0.12 

0.08±0.04 

0.53±0.09 

2.774=0.26 


23.454=0.45 



BX691 

17:01:05.996 

64:12:10.271 

2.190 

25.334=0.16 

0.224=0.10 

0.664=0.19 

4.644=0.17 

21.894=0.10 

21.654=0.10 

21.444=0.12 

22.024=0.27 

BX1118 

17:01:08.914 

64:16:12.554 

2.214 

24.554=0.16 

0.074=0.08 

0.574=0.11 

3.144=0.36 


23.674=0.29 



BX1047 

17:00:23.369 

64:15:31.444 

2.231 

24.554=0.16 

0.434=0.11 

1.044=0.16 

4.014=0.16 


22.024=0.10 


22.614=0.31 

BX1084 

17:01:32.349 

64:15:50.699 

2.245 

23.764=0.12 

0.454=0.07 

1.054=0.11 

3.894=0.10 


22.144=0.46 



BX931 

17:00:25.497 

64:14:29.805 

2.246 

25.154=0.16 

0.474=0.12 

1.124=0.28 

3.604=0.28 


22.834=0.11 



BX794 

17:00:47.301 

64:13:18.702 

2.253 

23.60±0.12 

0.354=0.06 

0.58±0.09 

3.104=0.18 


23.674=0.25 



BX1030 

17:00:46.895 

64:15:23.275 

2.285 

25.114=0.16 

0.474=0.12 

1.174=0.28 

4.024=0.22 


22.324=0.18 



BX709 

17:00:54.761 

64:12:19.973 

2.285 

25.184=0.16 

0.05±0.10 

0.77±0.15 

3.064=0.59 

23.37±0.38 

23.064=0.28 



BX898 

17:01:05.284 

64:14:06.570 

2.286 

24.484=0.14 

0.184=0.06 

0.574=0.11 

2.374=0.51 


22.174=0.21 



MD69 

17:00:47.610 

64:09:44.78 

2.289 

24.854=0.16 

0.374=0.09 

1.504=0.33 

4.874=0.12 

21.134=0.03 

20.894=0.04 

20.594=0.07 

20.634=0.11 

BX810 

17:01:31.117 

64:13:29.357 

2.290 

24.684=0.16 

0.364=0.09 

0.87±0.18 

3.994=0.19 

22.924=0.23 

22.784=0.23 



MD157 

17:00:52.191 

64:15:29.245 

2.293 

24.354=0.14 

0.354=0.07 

1.654=0.14 

4.464=0.10 


21.784=0.10 


21.574=0.12 

MD109 

17:01:04.482 

64:12:09.293 

2.294 

25.46±0.16 

0.26±0.10 

1.444=0.37 

3.604=0.28 

23.124=0.11 

22.864=0.10 



BX563 

17:01:15.755 

64:10:27.088 

2.296 

23.824=0.12 

0.424=0.07 

0.864=0.09 

3.154=0.18 


21.804=0.28 



BX939 

17:00:52.915 

64:14:36.023 

2.296 

24.464=0.14 

0.384=0.07 

0.694=0.13 

3.464=0.19 


23.054=0.13 



BX984 

17:01:02.164 

64:14:56.440 

2.299 

23.514=0.12 

0.174=0.04 

0.72±0.06 

3.444=0.14 


22.234=0.10 



BX919 

17:00:33.909 

64:14:21.626 

2.301 

24.434=0.14 

0.034=0.06 

0.654=0.10 

2.704=0.41 


23.514=0.23 



BX917 

17:01:16.110 

64:14:19.801 

2.302 

24.43±0.14 

0.28±0.07 

0.95±0.13 

4.044=0.12 


22.054=0.13 


22.194:0.37 

BX1133 

17:01:03.294 

64:16:19.717 

2.303 

24.384=0.14 

0.144=0.06 

0.82±0.13 

3.334=0.23 


22.854=0.10 



BX879 

17:01:17.768 

64:13:55.116 

2.308 

23.504=0.12 

0.614=0.04 

1.404=0.19 

2.844=0.22 

22.224=0.34 

22.224=0.15 



BX505 

17:00:48.224 

64:10:05.862 

2.309 

25.174=0.16 

0.454=0.12 

1.284=0.26 

4.364=0.19 

22.234=0.10 

22.114=0.10 

22.494=0.17 

22.104=0.15 

BX585 

17:01:13.491 

64:10:58.829 

2.311 

24.714=0.16 

0.124=0.08 

0.764=0.13 

3.154=0.36 



23.564=0.41 


MD103 

17:01:00.209 

64:11:55.576 

2.315 

24.234=0.14 

0.464=0.09 

1.494=0.27 

3.974=0.12 

21.324=0.10 

21.234=0.10 

21.344=0.14 


BX532 

17:00:40.154 

64:10:22.537 

2.331 

23.764=0.12 

0.224=0.06 

0.57±0.09 

3.634=0.12 

21.024=0.10 

21.104=0.10 

21.214=0.16 

21.564=0.21 

MD94 

17:00:42.020 

64:11:24.224 

2.333 

24.72±0.16 

0.944=0.15 

2.06±0.25 

5.154=0.15 

20.594=0.08 

20.434=0.03 

20.354=0.22 

20.284=0.14 

BX772 

17:00:53.076 

64:13:04.566 

2.346 

24.964=0.16 

0.20±0.09 

0.57±0.14 

3.294=0.31 


23.724=0.25 



MD174 

17:00:54.542 

64:16:24.760 

2.347 

24.564=0.16 

0.32±0.09 

1.504=0.27 

4.744=0.11 


20.514=0.10 


18.51 4=0.10 

BX490 

17:01:14.831 

64:09:51.691 

2.403 

22.884=0.06 

0.364=0.03 

0.92±0.05 

2.884=0.10 

22.054=0.10 

22.084=0.10 

21.954=0.11 

22.134=0.30 

BX581 

17:01:02.726 

64:10:51.299 

2.406 

23.874=0.12 

0.284=0.06 

0.62±0.10 

3.214=0.15 

23.114=0.15 

23.174=0.14 

22.604=0.21 

21.944=0.19 

BX759 

17:00:59.545 

64:12:55.455 

2.418 

24.434=0.14 

0.364=0.07 

1.294=0.18 

3.244=0.23 

21.814=0.31 

22.394=0.29 



BX561 

17:01:04.180 

64:10:43.834 

2.426 

24.654=0.16 

0.194=0.08 

1.044=0.16 

4.634=0.11 

21.544=0.10 

21.354=0.10 

21.194=0.10 

20.854=0.10 

BX575 

17:01:03.341 

64:10:50.932 

2.431 

23.824=0.12 

0.164=0.04 

0.74±0.06 

2.604=0.26 

23.104=0.12 

23.004=0.10 

23.034=0.41 










TABLE 2 


Continued 


Name 

a( J2000) 

<5( J2000) 

z 

TZab 

(G - 1Z ) AB 

(U n - G) AB 

n - K s b 

(IRAC-I)ab 

(IRAC-2)ab 

(IRAC-3)ab 

(IRAC-4)ab 

BX717 

17:00:56.995 

64:12:23.763 

2.438 

24.78±0.16 

0.34±0.09 

0.47±0.13 

2.73±0.54 


23.89±0.24 



BX826 

17:01:17.163 

64:13:41.960 

2.438 

24.45±0.14 

0.02±0.06 

0.52±0.09 

2.94±0.28 

22.65±0.25 

22.55=h0.10 



BX523 

17:00:41.712 

64:10:14.884 

2.471 

24.51±0.16 

0.46±0.11 

1.28±0.18 

3.52±0.25 

22.74±0.20 

22.68=h0.15 


22.07±0.35 

BX553 

17:01:02.200 

64:10:39.390 

2.474 

24.49±0.14 

0.14±0.06 

0.44±0.11 

2.77±0.41 

23.13±0.13 

23.18±0.19 



MD72 

17:00:42.777 

64:10:05.346 

2.533 

24.63±0.16 

0.45±0.11 

1.51±0.00 

3.42±0.25 

22.44±0.15 

22.30±0.20 



BX568 

17:01:28.789 

64:10:46.682 

2.537 

23.93±0.12 

0.38±0.06 

0.94±0.09 

3.20±0.15 

22.75±0.10 

22.81±0.10 

22.82±0.23 

22.87±0.25 

BX588 

17:01:01.469 

64:11:00.301 

2.539 

24.67±0.16 

0.53±0.11 

1.17±0.19 

3.22±0.31 

24.22±0.29 



22.97±0.29 

BX609 

17:01:9.988 

64:11:15.663 

2.573 

24.12±0.14 

0.56±0.09 

1.33±0.18 

3.06±0.27 

23.10±0.12 

22.94±0.13 

23.58±0.46 

23.16±0.38 

MD154 

17:01:38.387 

64:14:57.372 

2.629 

23.23ih0.ll 

0.73±0.04 

1.91±0.22 

3.55±0.09 

21.50±0.27 




BX591 

17:00:50.179 

64:11:01.162 

2.686 

24.68±0.16 

0.22±0.09 

1.06±0.16 

2.93±0.54 

23.48±0.14 

23.91±0.23 



MD128 

17:00:30.292 

64:13:15.165 

2.713 

23.41±0.11 

0.29±0.04 

1.24±0.08 

2.21±0.28 


23.67±0.16 



MD98 

17:01:04.387 

64:11:49.543 

2.738 

24.54±0.16 

0.70±0.10 

1.82±0.46 

2.94±0.54 

22.62±0.10 

22.94±0.22 



BX928 

17:01:19.623 

64:14:25.911 

2.755 

23.77±0.12 

0.66±0.07 

1.54±0.19 

2.95±0.22 


22.99±0.21 



MD152 

17:00:31.085 

64:14:54.048 

2.791 

23.43±0.11 

0.86±0.05 

1.92±0.29 

2.92±0.15 


21.62±0.10 


21.65=h0.11 

BX529 

17:00:44.618 

64:10:19.001 

2.875 

24.27±0.14 

0.25±0.06 

0.28±0.06 

2.74±0.41 

23.35±0.25 




MD155 

17:00:47.888 

64:15:05.446 

2.875 

24.03±0.14 

0.56±0.09 

1.83±0.43 

2.32±0.51 


23.31=h0.14 



MD126 

17:00:51.610 

64:13:05.346 

2.898 

23.94±0.12 

0.67±0.04 

2.01±0.46 

2.32±0.42 


22.82±0.15 




a Includes all galaxies with spectroscopic redshifts that are significantly detected in the ground-based A's-band image and in one or more of the IRAC bands. 
b 'IZ(AB) — K s (Ve ga), to make the comparison with other samples easier. This color can be converted to the AB system used in the fitting and in the plots in Figure l5l 
by subtracting 1.82. 


to 
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TABLE 3 
Modeling Results 


to 

00 


Object 

z 

log M*(CSF) a 

E(B-V) 

SFR b 

Age c 

log M*(r) d 

E(B-V) 

SFR e 

Age f 

T S 

BX627 

1.478 

9.42±0.36 

0.05 

7 

404 

9.42=1=0.43 

0.05 

7 

404 

const 

BX738 

1.710 

10.40±0.09 

0.10 

20 

1278 

10.35=1=0.05 

0.00 

4 

404 

100 

BX756 

1.738 

9.76i0.05 

0.34 

377 

15 

9.76=1=0.05 

0.34 

377 

15 

const 

BX681 

1.740 

10.50i0.24 

0.21 

246 

128 

10.65=1=0.03 

0.08 

1 

81 

10 

BX843 

1.771 

9.89±0.23 

0.21 

24 

321 

10.06=1=0.14 

0.04 

1 

255 

50 

BX470 

1.841 

10.814=0.09 

0.30 

64 

1015 

10.75=1=0.07 

0.23 

24 

321 

100 

BX846 

1.843 

10.444=0.10 

0.12 

27 

1015 

10.39=1=0.08 

0.03 

7 

360 

100 

BX574 

1.846 

10.164=0.19 

0.24 

51 

286 

10.34=1=0.10 

0.06 

3 

255 

50 

BX1043 

1.847 

9.914=0.25 

0.10 

32 

255 

9.94=1=0.10 

0.00 

0 

90 

10 

BX1087 

1.871 

10.244=0.11 

0.15 

68 

255 

10.31=1=0.08 

0.04 

3 

114 

20 

BX1032 

1.883 

10.084=0.16 

0.20 

53 

227 

10.12=1=0.11 

0.10 

4 

102 

20 

BX782 

1.930 

10.394=0.11 

0.10 

24 

1015 

10.33=1=0.07 

0.00 

6 

360 

100 

BX530 

1.942 

9.844=0.21 

0.34 

699 

10 

10.29=1=0.09 

0.08 

1 

72 

10 

BX807 

1.971 

9.864=0.16 

0.10 

10 

719 

9.85=1=0.14 

0.01 

3 

321 

100 

BX903 

1.976 

10.104=0.16 

0.14 

22 

571 

10.lli0.20 

0.10 

13 

360 

200 

BX536 

1.977 

10.654=0.06 

0.23 

138 

321 

10.74i0.06 

0.09 

20 

203 

50 

BX412 

1.981 

10.244=0.18 

0.20 

61 

286 

10.28i0.15 

0.12 

16 

161 

50 

BX526 

2.018 

10.064=0.15 

0.04 

16 

719 

9.97i0.18 

0.00 

8 

255 

100 

BX1075 

2.061 

8.984=0.31 

0.30 

115 

8 

8.97i0.31 

0.30 

88 

7 

10 

BX625 

2.078 

10.044=0.12 

0.07 

12 

905 

10.02i0.10 

0.01 

6 

454 

200 

BX592 

2.099 

9.304=0.22 

0.37 

199 

10 

9.82i0.24 

0.12 

2 

102 

20 

BX535 

2.113 

9.484=0.24 

0.43 

199 

15 

10.00i0.29 

0.23 

1 

72 

10 

BX1012 

2.115 

10.104:0.20 

0.07 

25 

509 

10.06i0.19 

0.01 

10 

255 

100 

BX691 

2.190 

11.044:0.03 

0.29 

40 

2750 

11.14i0.05 

0.13 

9 

2750 

1000 

BX1118 

2.214 

9.974=0.20 

0.07 

13 

719 

9.86i0.20 

0.00 

4 

286 

100 

BX1047 

2.231 

10.714=0.12 

0.26 

50 

1015 

10.67i0.07 

0.12 

8 

404 

100 

BX1084 

2.245 

11.094=0.11 

0.26 

120 

1015 

10.92i0.09 

0.10 

15 

404 

100 

BX931 

2.246 

10.194=0.27 

0.27 

38 

404 

10.24i0.27 

0.20 

15 

255 

100 

BX794 

2.253 

9.494=0.36 

0.29 

430 

7 

9.45i0.41 

0.28 

265 

7 

10 

BX1030 

2.285 

10.654=0.17 

0.27 

39 

1139 

10.58i0.17 

0.23 

22 

454 

200 

BX709 

2.285 

10.374=0.17 

0.10 

10 

2400 

10.37i0.27 

0.10 

10 

2400 

const 

BX898 

2.286 

10.734=0.11 

0.10 

20 

2750 

10.73i0.21 

0.10 

20 

2750 

const 

MD69 

2.289 

11.464=0.02 

0.35 

105 

2750 

11.51i0.05 

0.24 

40 

2200 

1000 

BX810 

2.290 

10.514=0.12 

0.20 

29 

1139 

10.41i0.09 

0.01 

3 

454 

100 

MD157 

2.293 

11.034:0.07 

0.28 

66 

1609 

10.89i0.07 

0.21 

22 

360 

100 

MD109 

2.294 

10.484:0.17 

0.17 

14 

2100 

10.48i0.25 

0.17 

14 

2100 

const 

BX563 

2.296 

10.574:0.23 

0.21 

92 

404 

10.57i0.28 

0.21 

92 

404 

const 

BX939 

2.296 

10.234=0.15 

0.16 

27 

640 

10.28i0.ll 

0.00 

3 

404 

100 

BX984 

2.299 

10.654=0.07 

0.12 

44 

1015 

10.55i0.05 

0.01 

10 

360 

100 

BX919 

2.301 

9.844=0.21 

0.07 

17 

404 

9.74i0.26 

0.00 

5 

161 

50 

BX917 

2.302 

10.964=0.07 

0.18 

33 

2750 

10.86i0.09 

0.14 

22 

1015 

500 

BX1133 

2.303 

10.374:0.12 

0.10 

20 

1139 

10.27i0.ll 

0.01 

5 

360 

100 

BX879 

2.308 

10.034:0.10 

0.42 

1473 

7 

10.02i0.13 

0.41 

1380 

7 

100 

BX505 

2.309 

10.794=0.13 

0.28 

38 

1609 

10.73i0.08 

0.20 

16 

571 

200 

BX585 

2.311 

10.034=0.28 

0.11 

17 

640 

9.99i0.26 

0.01 

4 

321 

100 

MD103 

2.315 

11.074=0.13 

0.28 

103 

1139 

10.96i0.13 

0.21 

38 

321 

100 

BX532 

2.331 

11.264=0.03 

0.14 

65 

2750 

11.33i0.08 

0.09 

36 

2750 

2000 

MD94 

2.333 

11.524=0.10 

0.49 

324 

1015 

11.43i0.10 

0.42 

113 

321 

100 

BX772 

2.346 

10.074=0.18 

0.10 

12 

1015 

10.06i0.16 

0.03 

5 

509 

200 

MD174 

2.347 

11.594:0.04 

0.34 

140 

2750 

11.65i0.05 

0.20 

40 

2500 

1000 

BX490 

2.403 

9.924=0.12 

0.29 

835 

10 

9.88i0.33 

0.28 

506 

9 

10 

BX581 

2.406 

10.374=0.11 

0.11 

33 

719 

10.36i0.10 

0.10 

26 

509 

500 

BX759 

2.418 

10.524=0.22 

0.23 

73 

454 

10.52i0.28 

0.23 

73 

454 

const 

BX561 

2.426 

11.314=0.02 

0.29 

79 

2600 

11.39i0.06 

0.16 

25 

2400 

1000 

BX575 

2.431 

9.954=0.20 

0.11 

44 

203 

9.99i0.25 

0.07 

25 

161 

100 



TABLE 3 


Continued 


Object 

z 

log M*(CSF) a 

E(B-V) 

SFR b 

Age c 

log M*(r) d 

E(B-V) 

SFR® 

Age f 

T S 

BX717 

2.438 

9.85=t0.38 

0.09 

14 

509 

9.99±0.42 

0.01 

6 

454 

200 

BX826 

2.438 

10.63±0.04 

0.04 

17 

2600 

10.69±0.06 

0.01 

11 

2300 

2000 

BX523 

2.471 

10.31±0.30 

0.25 

71 

286 

10.39±0.28 

0.18 

28 

227 

100 

BX553 

2.474 

10.34±0.10 

0.04 

14 

1609 

10.30±0.10 

0.01 

10 

806 

500 

MD72 

2.533 

10.61±0.28 

0.22 

56 

719 

10.61±0.30 

0.22 

56 

719 

const 

BX568 

2.537 

10.27±0.27 

0.16 

52 

360 

10.33±0.06 

0.01 

5 

227 

50 

BX588 

2.539 

9.47±0.22 

0.35 

295 

10 

9.44±0.27 

0.34 

210 

10 

20 

BX609 

2.573 

9.68±0.15 

0.35 

576 

8 

10.17±0.25 

0.00 

0 

128 

10 

MD154 

2.629 

10.94±0.24 

0.35 

677 

128 

10.96±0.09 

0.01 

6 

286 

50 

BX591 

2.686 

10.10±0.19 

0.06 

16 

806 

10.09±0.16 

0.00 

7 

454 

200 

MD128 

2.713 

9.38±0.09 

0.18 

287 

8 

9.37=t0.11 

0.18 

223 

7 

10 

MD98 

2.738 

9.82±0.26 

0.36 

436 

15 

10.47±0.27 

0.07 

6 

227 

50 

BX928 

2.755 

9.85±0.13 

0.35 

972 

7 

9.87±0.21 

0.36 

895 

6 

10 

MD152 

2.791 

10.28±0.06 

0.41 

2087 

9 

10.28±0.16 

0.41 

2087 

9 

const 

BX529 

2.875 

10.33±0.20 

0.03 

19 

1139 

10.28±0.24 

0.01 

15 

640 

500 

MD155 

2.875 

9.52±0.18 

0.26 

333 

10 

9.48=t0.25 

0.24 

202 

9 

10 

MD126 

2.898 

9.72±0.09 

0.29 

520 

10 

9.69=t0.23 

0.28 

380 

10 

20 


a Stellar mass and uncertainty, in solar units, for constant star formation model. 
b Star formation rate, in Mq yr —1 , for constant star formation model. 
c Age, in Myr, for constant star formation model. 

d Stellar mass and uncertainty, in solar units, for the best-fitting r model. The stellar mass uncertainty reflects the uncertainty in r. 
e Star formation rate, in M SU n yr — ■*■, for the best-fitting r model. 
f Age, in Myr, for the best-fitting r model. 

g Best-fitting star-formation decay timescale (SFR oc e —T ), in Myr. 


to 

CD 



TABLE 4 

Stellar Mass Estimates as a Function of K s Magnitude 


Sample 

(M*(r)) a 

(. M,{CSF)) b 

(M*(r)>(no IRAC) C 

All (72) 

10.32 ± 0.51 

10.30 ±0.53 

10.37 ±0.55 

K s < 21 (38) 

10.56 ±0.50 

10.55 ±0.51 

10.67 ±0.44 

K s < 20.6 (20) 

10.75 ±0.55 

10.75 ±0.56 

10.88 ±0.40 

Ks < 20 (9) 

10.96 ±0.55 

10.96 ±0.55 

11.01 ±0.40 


Mean and standard deviation of the log of the inferred stellar mass for the best-fit r models. 

Mean and standard deviation of the log of the inferred stellar mass for the CSF models. 

Mean and standard deviation of the log of the inferred stellar mass for the best-fit r model, excluding the IRAC photometry. 
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